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CHAPTER  I.  INTRODUCTION 

One  aspect  of  turbulent  shear  flows  that  needs 
empiricism  to  be  able  to  predict  it  is  the  turbulent  trans- 
port and  spreading  of  these  types  of  flows.  A turbulent 
plane  jet,  for  example,  has  a linear  widening  rate  away 
from  the  mouth  of  the  jet  which  has  not  been  predicted  from 
turbulent  flow  theories,  but  rather  determined  from  mea- 
surements themselves.  Various  models  have  been  proposed 
to  try  to  explain  the  spreading  rate  of  turbulent  develop- 
ing shear  flows.  In  recent  years,  the  recognition  of  some 
coherent  structures  in  turbulence  has  led  to  different 
models  and  approaches  that  attempt  to  explain  turbulent 
flows.  Some  of  these  were  aptly  summarized  by  Davies 
and  Yule  ( 1 974  )*  fol  1 owi  ng  a conference  in  Southampton. 

This  idea  of  coherent  structures  in  turbulent  flows 
is  an  extension  and  companion  to  the  concept  of  a distinct 
surface  separating  the  turbulent  primary  flow,  as  for 
instance  in  a jet,  from  the  entrained  external  nonturbu- 
lent  flow.  The  existence  of  this  interface  has  been  well 
defined,  and  will  be  alluded  to  in  future  paragraphs.  In 
the  last  two  or  three  years,  it  has  been  recognized  that 
occasionally  this  interface  separating  the  turbulent  from 
the  nonturbulent  flows  tends  to  fold-over,  or  seemingly 
Also  published  t JFM  (1975),  Vol.  69,  Part  3,  pp.  513-537. 
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engulf  the  external  nonturbulent  flow  in  a "nibbling" 
fashion.  In  the  presence  of  fold-over,  a region  of  the 
turbulent  fluid  may  have  nontu rbu 1 ent  f 1 u i d between  it  and 
the  main  fully  turbulent  region. 

The  overall  objective  of  this  work  was  to  try  to  gain 
further  understanding  into  the  structure  of  turbulent  flows, 
and  to  observe  the  effects  of  an  acoustic  disturbance  upon 
them.  The  specific  objective  of  this  work,  to  which  the 
attention  is  now  directed,  was  a measure  of  the  fold-over 
of  the  interface  between  the  turbulent  and  nonturbulent 
flows,  and  the  effect  of  an  acoustic  disturbance  upon  it. 

In  addition  to  that,  the  i ntermi ttency  and  crossing  fre- 
quency of  the  interface  were  recorded.  The  frequency  of 
the  acoustic  disturbance  was  purposely  chosen  to  match  that 
frequency  at  which  the  most  interaction  is  expected,  as 
exhibited  by  the  noticeable  increase  in  the  widening  rate 
of  the  acoustically  disturbed  jet  compared  with  the  jet 
without  any  sound  applied  to  it.  This  effect  of  the 
acoustic  field  on  the  jet  is  discussed  further  in  Chapter 
II. 

The  initial  goal  of  this  work  was  then  to  see  if  an 
acoustic  disturbance  in  the  form  of  a plane  traveling  wave 
at  specific  frequencies  could  indeed  change  the  occurrence 
of  fold-over  in  a two  dimensional  plane  turbulent  jet. 

These  specific  frequencies  will  be  discussed  further  in 
Chapter  II.  Folding  of  the  interface  as  discussed  by 
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Townsend  (1970),  Kohan  (1969),  Paizis  (1972),  LaRue  (1973), 
LaRue  and  Libby  (1974),  Davies  and  Yule  (1974),  Paizis  and 
Schwarz  (1974),  and  Mulej  (1975)  is  defined  as  the  process 
by  which  a region  of  turbulent  fluid  has  nonturbulent 
fluid  between  it  and  the  fully  turbulent  region.  In  an 
attempt  to  describe  this  folding  phenomenon,  some  investi- 
gators have  used  different  terms.  LaRue  (1973)  for 
example,  used  "overhangs,"  Mulej  (1975)  used  "folding" 
and  "fold-over"  interchangeably,  and  Paizis  and  Schwarz 
(1974)  referred  to  it  as  simply  "folding."  This  will  be 
discussed  further  in  Chapter  III.  The  term  "fold-over," 
will  be  used  in  the  present  work. 

Mulej  ( 1 975  ) measured  the  percent  fold-over,  4'^  in  a 
plane  turbulent  jet,  and  found  that  it  decreased  with  the 
lateral  position,  y/b,  of  the  jet.  (The  description  of 
the  jet  is  given  in  Chapter  III.)  The  present  research 
was  partly  motivated  by  his  report.  A new  and  hopefully 
better  method  of  measuring  fold-over  was  sought.  The 
method  used  by  Mulej  (1975)  is  discussed  and  evaluated  in 
Chapter  III.  The  method  adopted  in  this  research  is  a 
digital  method.  Analog  signals  were  digitized  and  processed 
in  a digital  computer. 

Fold-over  measurements  required  detecting  the  turbu- 
1 ent-nonturbu 1 ent  interface,  using  two  hot  wire  probes 
separated  by  a pre-selected  distance  in  the  lateral 
direction  of  the  flow.  Corrsin  (1943)  was  credited  with 
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discovering  turbulent  shear  flows  may  exhibit  flows 

which  are  intermittently  turbulent  or  nonturbu 1 ent . He 
was  exploring  the  flow  field  of  a heated  jet.  The  exist- 
ence of  this  phenomenon  was  later  established  for  all  free 
shear  flows,  and  Townsend  (1948)  proposed  a measure  of  the 
degree  of  i ntermi ttency . This  became  known  as  the  inter- 
mittency  factor,  y,  and  is  defined  as  the  fraction  of  time 
the  flow  is  turbulent.  Tne  i ntermi ttency  function,  I(t), 
was  also  due  to  Townsend  (1949).  This  function  is  a ran- 
dom square  wave  which  is  equal  to  1 (unity),  when  the  flow 
is  turbulent,  and  zero  when  it  is  nonturbu 1 ent . 

Recognizing  the  difficulty  involved  in  detecting  the 
tu rbu 1 en t-non turbu  1 en t interface,  the  objectives  of  this 
research  were  expanded  and  are  three-fold:  (1)  to  detect 
the  t u rbu 1 ent -non tu rbu 1 en t interface  using  both  analog 
and  digital  techniques,  (2)  to  measure  fold-over  using  the 
developed  detection  technique,  and  (3)  to  observe  the  ef- 
fects of  a plane  traveling  acoustic  wave  of  selected  fre- 
quencies on  fold-over,  crossing  frequency  of  the  interface, 
and  i ntermi ttency.  The  interface  crossing  frequency  is 
defined  as  the  rate  at  which  the  i ntermi ttency  function, 
I(t),  changes  from  unity  to  zero  or  vice  versa. 

Whereas  some  investigators  contend  that  fold-over  is 
an  important  part  of  the  entrainment  process,  only  very 
few  attempts  have  been  made  to  quantify  it.  In  the  plane 
turbulent  jet,  for  example,  the  only  reported  attempt  to 


5 


measure  fold-over  was  made  by  Mulej  (1975).  Some  of  the 
questions  raised  by  his  work  are:  (1)  Why  does  the  de- 
tection of  fold-over  seem  to  increase  as  one  moves  towards 
the  center-line  of  the  jet  where  the  flow  is  fully  turbu- 
lent? (2)  Do  fold-over  events  actually  increase  as  one 
moves  towards  the  center-line  of  the  jet?  (3)  Is  measur- 
ing the  percent  fold-over  the  best  way  to  quantify  fold- 
over? 

This  work  is  then  a natural  extension  of  the  work  done 
by  Jenkins  (1974)  and  Mulej  (1975)  in  plane  turbulent  jets 
without  the  application  of  an  acoustic  disturbance.  The 
effects  of  an  acoustic  wave  on  turbulent  jets  have  been 
reported  for  many  years.  None  of  the  previous  works  con- 
sidered the  intermittent  region  of  the  jet.  Simcox  (1969), 
for  examp  1 e,  pred i c ted  theoretically  that  the  spreading 
rate  of  the  jet  is  enhanced  by  an  acoustic  disturbance. 

It  was  proposed  by  him  that  this  interaction  of  an  acoustic 
wave  with  the  jet  may  be  applicable  to  the  development  of 
a new  type  of  amplifier.  Experiments  by  Kaiser  (1971)  and 
others  confirmed  that  the  widening  rate  of  the  jet  was 
increased  by  the  application  of  an  acoustic  wave.  Research 
to  establish  the  mechanism  through  which  an  applied  sound 
interacts  with  a plane  turbulent  jet  is  being  conducted 
by  Chambers  (1976).  The  location  in  the  jet  at  which  this 
interaction  occurs  is  also  being  investigated  by  him. 
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It  is  hoped  that  the  present  research  will  provide 
more  information  concerning  the  interaction  of  an  acoustic 
wave  with  a plane  turbulent  jet,  as  well  as  the  structure 
of  the  t ur bu 1 en t - non tu r bu 1 ent  interface.  This  might  lead  ^ 

to  a better  understanding  of  turbulence  and  its  interaction 
with  an  acoustic  disturbance.  i 

i 
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CHAPTER  II.  THE  TWO-DIMENSIONAL  JET 

1 . Characteristics  of  the  Jet 
Among  the  simplest  types  of  turbulent  shear  flows  is 
a jet.  It  is  an  unbounded  flow  and  is  formed  when  a fluid 
moving  with  a velocity  greater  than  that  of  the  surrounding 
medium  discharges  through  an  orifice  into  the  medium.  (The 
medium  considered  here  is  air  and  the  fluid  forming  the 
jet  is  also  air.)  The  resulting  diffusion  pattern  can  be 
divided  into  three  regions:  the  initial  region,  the 
transitional  region,  and  the  main  region  (see  Figure  2-1). 

Within  the  initial  region  of  the  jet,  there  exists  a 
potential  core  bounded  by  two  mixing  layers.  The  velocity 
in  this  potential  core  remains  constant.  At  the  mouth  of 
the  jet,  this  velocity  is  equal  to  the  discharge  velocity. 

As  the  flow  spreads,  the  axial  velocity  decreases. 
Turbulent  flow  results  due  to  the  gradient  and  corresponding 
exchange  caused  by  the  initial  velocity  discontinuity  be- 
tween the  jet  and  the  surrounding  media.  The  ensuing 
turbulent  stresses  bring  about  a transverse  exchange  of 
momentum  between  the  fluid  and  the  surrounding  media.  The 
"eating  up"  of  the  jet  beyond  the  initial  area  shows  up 
both  in  its  widening  as  well  as  in  variation  of  the  velocity 
along  its  axis  (see  Abramovich  (1963). 
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The  mean  velocity  profile  in  the  transitional  region 
will  essentially  be  fully  developed;  however,  a finite 
distance  is  required  for  the  turbulence  to  become  fully 
developed. 

The  main  region  begins  when  the  turbulence  becomes 
fully  developed.  It  is  c ha r ac t er i zed  by  a continuous  ac- 
celeration of  the  surrounding  fluid  accompanied  by  a 
deceleration  of  the  fluid  in  the  central  region.  The  mea- 
surements presented  later  were  taken  in  this  region  of  the 
flow  field. 

There  exists  a zone  of  demarcation  between  turbulent 
and  nonturbulent  regions  of  the  shear  flow.  This  zone  is 
characterized  by  a sharp,  randomly  convoluted  interface. 

The  turbulent-nonturbulent  region  (see  Figure  2-2)  is  the 
intermittent  region  with  which  this  work  is  concerned. 

The  shape  of  the  jet  profile  remains  unchanged  when 
an  acoustic  disturbance  is  applied  to  it.  However,  the 
spreading  rate  of  the  jet  is  increased  and  the  geometric 
virtual  origin  shifted  (see  Figure  2-3)  when  the  disturbance 
is  at  a specific  frequency. 

2.  General  Review  of  Previous  Work  Related  to 
turbulent  Jets 

It  was  observed  earlier  that  the  only  reported  attempt 
to  quantify  fold-over  in  a plane  turbulent  jet  was  that 
due  to  Mulej  (1975).  It  was  also  noted  that  no  work  has 
been  done  to  establish  the  effects  of  acoustic  disturbances. 


g Due  to  Acoustic  Disturbance, 
ser  (1971 )). 
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I at  specific  frequencies,  on  the  structure  of  the  turbulent- 

nonturbulent  interface  of  a two  d iniens  i oi\a  1 plane  turbulent 
jet.  This  section  therefore  aims  at  presenting  a brief 
review  of  the  literature  of  some  of  the  work  that  has 
been  done  by  previous  investigators  studying  turbulent  jets. 
Studies  in  other  flow  fields  like  wakes  and  boundary  layers 
will  be  mentioned  where  they  help  to  explain  concepts  like 
i nterm i t tency  and  conditional  sampling.  A specific  review 
of  the  previous  techniques  used  by  other  investigators  to 
detect  the  tu r bu 1 on t -non t u r bu 1 en t interface  will  be  given 
in  Chapter  111  when  the  detection  of  the  turbulent- 
nonturbulent  interface  is  discussed.  Also  the  work  done 
in  wakes  and  boundary  layers  to  measure  fold-over  will  be 
; ( presented  in  Chapter  111.  Hence  this  review  will  bo  di- 

vided into  two  parts:  (a)  the  undisturbed  flow  and  (b)  the 
; acoustically  excited  flow. 

t 

I 

a.  The  Undisturbed  Flow 

It  was  noted  in  Chapter  1 that  the  c ha-rac  t or  i s t i ca  1 1 y 
intermittent  hot  wire  signal  obtained  from  probes  placed  in 
the  t ur bu 1 en t -non t u r bu 1 en t region  of  a heated  jet  flow 
was  first  observed  by  Corrsin  (19A3).  The  existence  of 
the  phenomenon  of  i n t erm i t t ency  was  later  established  for 
all  free  shear  flows,  and  Townsend  (194H)  proposed  a mea- 
sure of  the  degree  of  i nterm  i t tency . This  became  known  as 
I the  i nt ermi t tency  factor,  >,  and  is  defined  as  the  fraction 

of  time  the  flow  is  turbulent.  The  1 ntermi t t ency  function. 


I(t),  was  also  due  to  Townsend  (1949).  This  function  is 
a random  square  wave  which  is  equal  to  1 (unity)  when  the 
probe  is  in  turbulent  flow  and  zero  when  it  is  in  the  non- 
turbulent  flow. 

Townsend  (1948)  first  explored  statistical  properties 
of  the  intermittent  region  for  the  two-dimensional  wake 
and  Corrsin  and  Kistler  (1954)  investigated  the  boundary 
layer.  The  thickness  of  the  interface  (or  superlayer,  as 
termed  by  Corrsin  and  Kistler  (1954))  was  investigated  by 
Corrsin  and  Kistler  (1954),  and  was  found  to  be  proportional 
to  the  Kolmogoroff  microscale  (v  /Cq)“.  The  existence 
of  this  interface  has  been  well  defined  following  measure- 
ments conducted  by  investigators  like  Kovasznay  (1967), 
Uberoi  and  Freymuth  (1969),  LaRue  (1973)  and  Kibens  (1968). 
Kovasznay  (1967),  argued  that  the  thickness  of  the  turbulent 
nonturbulcnt  interface  (superlayer)  is  of  the  same  order 
as  the  smallest  turbulent  eddies  and  estimated  this  thick- 
ness to  be  lOvVjj  ^ where  v is  the  kinematic  viscosity  and 
Vq  is  the  average  entrainment  velocity.  LaRue  (1973),  and 
also  LaRue  and  Libby  ( 1 9 74  ) argued  that  in  a wake  the  order  of  a 
minimum  duration  of  a fully  turbulent  "patch"  of  fluid  was 
equal  to  twice  the  Kolmogoroff  length  transit  time  or 
approximately  where  is  Kolmogoroff  length  and 

U is  the  free  stream  velocity.  Using  U obtained  by  Uberoi 
and  Freymuth  (1969)  they  obtained  the  minimum  duration  to 
be  160  VIS.  Further  discussion  of  the  interface  and  its 
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thickness  will  be  given  in  Chapter  III, 

Conventional  mean  measurement  techniques  require  con- 
tinuous averaging  of  flow  properties  as  obtained  from  a 
stationary  hot-wire  probe.  These  techniques  are  inherently 
incapable  of  examining  separately  the  detailed  properties 
of  the  turbulent  and  nonturbulent  flow,  or  of  the  trans- 
ition (or  interface)  between  the  two  states.  Corrsin  and 
Kistler  (1954)  realized  the  need  for  a new  technique.  The 
method  of  conditional  sampling  was  developed  primarily  by 
Kibens  and  Kovasznay  (1969)  in  answer  to  this  need.  Among 
the  early  workers  in  this  field  were  Kibens  (1968), 
Kovasznay,  Kibens  and  Blackwelder  (1970),  Jenkins  (1974) 
and  others.  Jenkins  (1974)  measured  the  velocity  and 
temperature  profiles  with  respect  to  the  interface  in  a 
heated  two-dimensional  plane  jet,  and  presented  an  up-to- 
date  literature  review  on  conditional  sampling. 

Phillips  (1972)  looked  at  the  part  played  by  the  inter- 
face in  the  entrainment  of  the  surrounding  ambient  fluid. 

A pseudo- Lag ra ng i a n description  of  the  entrainment  process 
was  proposed  by  him.  Paizis  and  Schwarz  (1974)  studied 
the  shape  of  the  interface  and  noted  the  significant  amount 
of  fold-over  which  occurred  in  the  interface.  They  also 
measured  the  characteristic  scales  and  convection  velocity 
of  the  interface  for  a turbulent  wall  jet.  More  discussion 
on  the  significant  amount  of  fold-over  measured  by  them 
is  made  in  Chapter  III. 


d 
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Mulej  (1975)  measured  the  lateral  component  of  the 
velocity  of  the  turbul ent-nonturbul ent  interface.  Pre 
liminary  measurements  of  the  amount  of  fold-over  in  the 
interface  as  a function  of  location  were  also  taken  by 
him. 


b.  The  Acoustically  Excited  Flow 

A survey  of  the  interaction  of  a transverse  acoustic 
wave  with  various  flows  is  presented  here.  This  discussion 
is  not  limited  to  plane  traveling  waves  and  is  not  intended 
to  be  an  exhaustive  one.  The  emphasis  is  on  turbulent  jet 
flows. 

Two  types  of  interaction  should  be  distinguished: 

(1)  that  interaction  that  promotes  an  earlier  transition 
in  the  case  of  an  initially  laminar  jet;  and  (2)  that  which 
causes  an  apparent  enhancement  of  mixing  in  the  fully 
turbulent  jet.  These  interactions  will  be  discussed  fur- 
ther in  this  chapter. 

The  influence  of  sound  on  a flow  field  has  been  known 
and  documented  for  many  years.  Chanaud  and  Powell  (1962) 
claim  that  the  first  reported  observation  of  this  phenome- 
non was  made  in  1838  by  LeConte,  who  noticed  that  gas  lights 
had  a tendency  to  flicker  when  in  the  presence  of  music. 

It  was  not  until  1 935  that  B^^own  ( 1 935  ) investigated  this 
phenomenon.  He  studied  sound  sensitive  laminar  jets,  and 
concluded  that  the  influence  of  sound  was  basically  an 
effect  upon  the  laminar-turbulent  transition. 


16 


Simcox  (1969)  first  represented  the  vorticity  fluctua- 
tions due  to  turbulence  by  a single  frequency  harmonic 
function,  and  superimposed  to  it  an  acoustic  disturbance. 

He  noticed  that  the  disturbance  caused  an  increase  in  the 
vorticity.  The  single  frequency  representation  was  then 
converted  to  a spectral  representation  using  Fourier  trans- 
forms which  described  a two-dimensional  turbulent  jet  flow. 
He  analyzed  the  interaction  of  a traveling  acoustic  wave 
with  the  turbulent  flow  field,  and  claimed  that  there  were, 
in  general,  two  types  of  interaction:  a direct  interaction 
(described  by  the  single  eddy  theory)  and  a secondary 
interaction  among  the  eddies  due  to  the  addition  of  energy 
by  the  acoustics.  His  analysis  suggested  an  increase  of 
the  turbulence  intensity,  and  no  change  of  the  Reynolds 
stress.  He  predicted  an  increase  in  the  turbulent  mixing, 
and  further  inferred  an  increased  spreading  rate  of  the  jet, 
since  the  spreading  was  proportional  to  the  kinetic  energy 
of  the  turbulence,  which  in  turn,  is  a function  of  the  tur- 
bulence intensity.  Although  the  theory  that  the  spreading 
rate  of  the  jet  is  proportional  to  the  kinetic  energy  of 
the  turbulence  has  not  been  proved  experimentally,  it  ap- 
pears to  be  the  consensus  among  other  investigators  like 
Abramovich  (1963)  and  Townsend  (1956).  Simcox  (1969)  fur- 
ther claimed  that  the  change  in  the  spreading  rate  of  the 
jet  can  be  used  to  produce  a control  on  the  switching  or 
proportioning  of  flow  in  a fluidic  device,  and  that  the 
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interaction  may  be  applicable  to  the  development  of  a new 
type  of  amplifier.  Part  of  the  work  by  Simcox  (1969)  is 
published  in  Simcox  and  Hoglund  (1971). 

Williams  (1969)  considered  a longitudinally  resonant 
acoustic  field  in  a developing  turbulent  pipe  flow,  and 
concluded  that  the  input  in  some  cases  led  to  a decrease  in 

the  microscales  and  macroscales  of  turbulence.  Although 
he  did  not  consider  jet  flows,  he  noticed  some  structural 
changes  in  turbulence  due  to  an  acoustic  signal. 

In  his  experimental  program  concerning  the  effect  of 
sound  on  a reattaching  jet  at  low  Reynolds  numbers,  Weinger 
(1965)  noticed  that  the  reattachment  point  moved  upstream 
with  increasing  acoustic  intensity  at  selective  frequencies. 
He  concluded  that  the  1 am i na r- tu r bu 1 ent  transition  shifted 
upstream  with  increasing  acoustic  intensity  and  that  once 
the  transition  point  reached  the  orifice,  the  reattachment 
point  began  to  move  downstream.  This  agreed  with  Brown's 
(1935)  observation. 

Roffman  and  Toda  (1969),  experimenting  with  both 
laminar  and  turbulent  jets,  found  an  increase  in  the  jet 
widening  rate  at  selective  frequencies.  They  also  found 
that  the  Strouhal  number  (St  = FD/Uq)  remained  constant 

over  the  range  of  their  investigation  (200  £ R-  9000). 

0 

It  was  Powell  (1953)  who  first  recognized  the  significance 
of  the  Strouhal  number  and  Reynolds  number  in  analyzing  the 
sensitivity  of  jets  to  periodic  disturbances. 


Glass  (1968)  undertook  an  experimental  investigation 
concerning  the  acoustic  feedback  on  supersonic  jets.  He 
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I 

i 

I 

placed  a plate  near  the  flow  field  in  order  to  reflect 

I 

the  aerodynami ca 1 1 y generated  sound  back  into  the  jet.  He  i 

found  that  this  type  of  acoustic  input  destroyed  the  shock 
structure  of  the  jet  and  caused  an  increase  of  spreading  j 

and  decay  rates.  The  effect  appeared  to  be  at  a maximum  J 

i 

at  10,000  Hz,  although  it  was  noticeable  at  other  j 

frequencies.  | 

Rockwell  (1971)  collected  most  of  the  data  of  the  j 

investigators  mentioned  above  and  presented  them  in  terms  | 

1 

of  two-dimensionless  groups:  the  Strouhal  number  and  the  j 

Reynolds  number.  His  plot  (Rockwell  (1971),  Figure  3)  j 

of  Strouhal  number  versus  Reynolds  number  on  a log-log  ‘j 

scale  showed  that  the  log  of  the  dimensionless  sensitive  ■ 

I 

frequencies  for  two-dimensional  jets  varies  approximately 

linearly  with  the  logarithm  of  the  Reynolds  number.  Hence,  ; 

it  is  possible  to  determine  the  frequency  at  which  a given 

flow  is  most  sensitive  to  an  acoustic  disturbance. 

The  influence  of  sound  on  turbulent  jet  flows 
may  be  through  an  effect  on  transition,  acoustic 
streaming,  an  induced  oscillation  due  to  aerodynamic 
loading  or  a direction  interaction  with  the  turbulent 
field  itself.  Simcox  (1969)  argued  qualitatively,  that 
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for  subsonic  air  jets  of  moderate  Reynolds  number,  neither 
acoustic  streaming,  nor  aerodynamic  loading  can  lead  to  a 
change  in  the  spreading  rate.  Experiments  conducted  by 

Kaiser  (1971)  confirm  Simcox's  (1969)  claim.  Hence 
the  noticeable  effects  of  an  acoustic  input  on  a free  jet 
are  through  an  upstream  shift  or  modulation  of  the  laminar- 
turbulent  transition,  and/or  a direct  interaction  with  the 
turbulent  flow  itself. 

To  evaluate  Simcox's  theoretical  predictions,  Kaiser 
(1971)  conducted  his  experiments  using  a two-dimensional 
turbulent  jet.  He  noted  that  the  similarity  of  the  flow 
remained  unaltered  when  an  acoustic  signal  was  introduced 
normal  to  the  jet.  Hence,  the  form  (U/0^^  = E(n))  remained 
unaltered.  He  confirms  that  the  jet  spreading  rate  is 
enhanced  by  an  acoustic  disturbance.  From  the  plot  of 
widening  rate  versus  frequency,  Kaiser  (1971)  noticed  a 
periodic  behavior  which  he  could  not  explain.  He  noticed 
that  as  the  Reynolds  number  of  the  flow  was  decreased,  the 
frequency  for  which  the  maximum  widening  rate  occurred 
also  decreased.  This  agreed  with  Rockwell's  (1971)  find- 
ings. The  geometric  origin  of  the  jet  tends  to  move  down- 
stream at  the  frequencies  where  the  maximum  widening  rate 
occurred.  The  velocity  origin  remained  approximately  zero. 

Some  of  the  issues  not  answered  by  Kaiser  (1971)  and 


i 


other  experimenters  were:  the  cause  of  the  periodic 
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behavior  of  the  widening  rate  as  a function  of  frequency, 
and  the  mechanism  through  which  an  applied  sound  field 
interacts  with  a plane  turbulent  jet.  Kaiser's  (1971) 
results  are  also  published  in  Goldschmidt  and  Kaiser  (1971). 

More  recent  studies  have  been  made  to  better  under- 
stand the  behavior  of  a disturbed  jet.  Chan  (1974)  studied 
spatial  waves  in  a low  speed  axisymmetric  turbulent  free 
jet.  His  results  showed  that  the  wave-numbers  of  the  pres- 
sure waves  increased  monotoni ca 1 ly , while  the  phase  veloci- 
ties decreased  with  the  jet  Strouhal  number.  Further 
studies  on  phase  velocities  were  made  by  Bechert  and 
Pfizenmaier  (1975).  They  applied  axial  disturbances  to  a 
free  jet  and  reported  the  presence  of  "ultra-fast"  phase 
velocities  which  exceeded  the  mean  velocity  of  the  jet. 
Hussain  and  Zaman  (1975)  studied  the  effect  of  an  acoustic 
excitation  on  the  turbulence  structure  of  a circular  jet 
operated  at  the  Reynolds  number  range  11,000-113,000.  The 
acoustic  excitation  was  in  a plenum  chamber.  They  observed 
that  at  a Strouhal  number  of  1.6,  the  turbulent  intensity 
in  the  near  field  of  the  nozzle  was  strongly  suppressed  but 
that  at  higher  Strouhal  numbers  the  effect  of  excitation 
was  insignificant.  Stiffler  (1975)  found  that  the  signal 
transit  velocities  in  a turbulent  plane  jet  were  much  less 
than  the  jet  center-line  velocity. 

Merkine  and  Liu  (1975)  studied  the  noise-producing 
large-scale  wavelike  eddies  in  a plane  turbulent  jet.  They 
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found  that  for  a given  frequency  parameter,  different  modes 
had  different  streamwise  lifetimes. 

It  was  pointed  out  earlier  that  this  survey  was  not  [ 

I 

intended  to  give  a historical  account  of  all  the  research  i 

studies  that  have  been  done  in  the  interaction  of  an  acous- 
tic wave  with  turbulent  jets.  It  aimed  at  pointing  out 
some  of  the  observations  that  have  been  made.  Whereas  many 
different  studies  have  been  conducted  concerning  the  inter- 
action of  an  acoustic  wave  with  turbulent  jets,  no  report 
has  been  made  of  the  effect  of  an  acoustic  disturbance  on 
i ntermi ttency , fold-over  and  crossing  frequency  of  the 
turbulent-nonturbulent  interface. 
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CHAPTER  III.  INVESTIGATION  TECHNIQUES  AND  EQUIPMENT 

1 . Introduction 

It  was  previously  stated  that  the  conventional  mean 
measurement  techniques  require  continuous  averaging  of 
flow  properties  as  obtained  from  a stationary  hot-wire 
probe.  Those  measurement  techniques  were  inherently  in- 
capable of  examining  separately  the  detailed  properties  of 
the  turbulent  and  nonturbulent  flows,  or  of  the  transition 
(or  interface)  between  the  two  states.  The  method  of 
conditional  sampling  was  developed  primarily  by  Kibens  and 
Kovasznay  (1969)  in  answer  to  this  problem.  Two  types  of 
conditional  averages  can  be  defined  within  the  turbulent 
region  of  a jet.  These  are  the  "zone"  averages  (taken  at 
a fixed  point  in  space  throughout  the  turbulent  history  of 
the  flow  at  that  point)  and  the  "point"  averages,  which  are 
taken  at  a point  within  the  turbulent  flow  and  averaged 
whenever  the  interface  is  crossing  some  other  reference 
point.  The  former  ("zone"  averages)  were  employed  in  this 
work.  The  continuous  signals  originating  from  two  hot-wire 
probes  placed  in  the  intermittent  zone  of  the  free  turbu- 
lent jet  were  classified  into  "turbulent"  or  "nonturbulent" 
by  digital  methods  and  processed  in  a computer.  Before  the 
details  of  this  method  are  presented,  a brief  discussion 
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of  the  experimental  setup  and  the  flow  field  will  be  given, 

2 . Experimental  Setup 

a.  The  Flow  Field 

The  two  dimensional  plane  jet  used  for  this  research 
was  essentially  the  same  as  that  used  by  Jenkins  (1974)  and 
Mulej  (1975)  (see  Figure  3-1).  It  was  powered  by  a 0.56  kW 
(0.75  hp)  blower.  A 4 kW  mesh  wire  heating  element  located 
in  a 56  X 51  cm  plenum  chamber  was  capable  of  heating  the 
flow.  This  heating  capability  used  by  Jenkins  (1974)  was 
not  used  for  this  work.  The  flow  was  directed  from  the 
plenum  chamber  via  a gradual  contraction  to  a smaller  5.98 
in  X 11.97  in. (15.2  x 30.4  cm)  plenum  chamber  having  flow 
straightening  elements.  It  was  then  discharged  through  a 
vertical  slot  0.5  i n.  x 11.97  in.  (1  . 27  x 30.4  cm)  on  a 35.98 
i n.  X 11.97  in.  (91.4  x 30.4  cm)  wall.  Two  confining  hori- 
zontal walls  35.98  i n.  x 47.64  in.  (91.4  x 121.0  cm)  were  built 
to  maintain  the  two  dimensionality  of  the  flow.  The 
Reynolds  number  (based  on  the  jet  slot  width)  was  held 
constant  at  1.78  x 10^. 

A two-dimensional  traversing  mechanism  was  used  to  move 
the  probe  appratus  (see  Figure  3-1).  The  range  of  the  tra- 
versing device  was  60  slot-widths  longitudinally  ( x -d i r ec t i on  ) 
and  44  slot-widths  laterally  (y-d i rec t i on  ) with  an  accuracy 
better  than  -0.02  in. (-0.04  cm).  An  L.C.  Smith  actuator 


i 

j 


that  permitted  independent  movement  of  one  of  the  sampling 
probes  was  mounted  on  the  probe  actuator  mount 
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(Figure  3-1).  The  other  sampling  probe  was  mounted  to,  and 
was  moved  manually,  with  the  mount.  The  actuator  provided 
electrical,  remote  controlled,  traversing  of  one  of  the 
sampling  probes.  A DC  voltage  across  a 10-turn  potentiometer 
provided  a remote  voltage  proportional  to  the  lateral  tra- 
verse. The  total  available  traverse  was  approximately  30 
slot- widths  relative  to  the  actuator  mount  with  an  accuracy 
of  -0.01  in. (-0.02  cm).  The  use  of  the  actuator  was  later 
abandoned  since  the  calibration  curve  (voltage  versus  the 
lateral  traverse)  varied  appreciably  before  and  after  a run. 
The  probe  separation  (in  the  lateral  direction)  was  now 
achieved  manually  and  this  was  as  accurate  as  reading  the 
scales  on  a ruler.  Only  the  top  probe  was  moved  manually 
(Figure  3-2).  The  probes  could  bo  moved  up  to  a separation 
of  8 slot-widths,  or  4 in,  (10.16  cm).  Most  measurements 
were  taken  at  a probe  separation  of  0.125  in,  (0.3175  cm). 
Figure  3-2  shows  the  two  hot  wire  probes.  The  hot  wire 
probes,  part  of  the  traversing  mechanism  and  the  confining 
horizontal  walls  of  the  jet  can  be  seen  in  Figure  3-3. 

The  block  diagrams  of  the  experimental  setup  used  to 
take  the  various  measurements  and  to  digitize  the  recorded 
analog  signals  are  shown  in  Figures  3-4  and  3-5,  respec- 
tively. Figure  3-6  shows  part  of  the  instrumentation. 

The  outputs  of  the  hot  wire  probes,  which  were  located  at 
preselected  stations  in  the  flow  field,  were  recorded 


Sample  Rate 
Generator 
Model  1572 


FM  Tape  Recorder 
Model  PI6104 

Figure  3-5.  Schematic  of  Digitization  Process 
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simultaneously.  Figure  3-7  shows  a typical  record  from  the 
visicorder.  The  time  derivatives  of  the  hot  wire  voltages 
were  recorded  with  an  FM  tape  recorder.  Figure  3-8  shows 
a typical  visicorder  output  of  the  derivative  signals. 
Comparing  this  figure  with  Figure  3-7,  it  can  be  seen  that 
the  derivative  signals  discriminate  between  the  turbulent- 
nonturbulent  zones  much  better  than  the  voltage  signals. 
The  choice  of  the  discrimination  signal  will  be  discussed 
in  detail  in  a later  section.  The  procedure  used  to  digi- 
tize the  time  derivatives  of  the  hot  wire  voltages  is  pre- 
sented in  Appendix  A. 

b.  Acoustic  Excitation 

The  acoustic  disturbance  was  generated  by  an  Altec 
290E  acoustic  driver,  connected  to  a rectangular  section 
exponential  horn,  designed  with  a low  frequency  cutoff  of 
166  Hz  (see  Figure  3-9).  The  exponential  horn  was  designed 
by  Chambers  (1  975).  The  horn  was  28  in.  (71.12  cm)  long 
and  had  a mouth  7.25  in- (18.415  cm)  by  19  in. (48.26  cm). 

The  input  to  the  driver  was  supplied  by  a Brliel  and  Kjaer 
Beat  Frequency  Oscillator  type  1022  coupled  to  a Marantz 
Model  9 Amplifier.  A Systron-Donner  Model  114  Counter 
measured  the  frequency  of  the  signal. 

It  was  noted  in  Chapter  II  that  there  is  a frequency 
at  which  a given  jet  is  most  sensitive  to  an  acoustic 
disturbance.  It  was  also  noted  there  that  Rockwell  (1971) 
collected  data  from  different  investigators  and  presented 


Figure  3>7.  Typical  Visicorder  Output  of  Hot  Wire  Voltages 
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them  in  terms  of  two  dimensionless  groups:  the  Strouhal 
number  and  the  Reynolds  number.  He  plotted  the  log  of  the 
Strouhal  number  versus  the  log  of  the  Reynolds  number. 

This  curve  was  approximately  a straight  line.  From  this 
curve,  it  is  possible  to  determine  the  frequency  at  which 
a given  flow  will  be  most  sensitive  to  an  acoustic  dis- 
turbance. Based  on  that  curve,  the  Strouhal  number  at 
which  the  jet  should  be  most  sensitive  to  the  acoustic 
disturbance  was  computed.  From  this,  the  first  frequency 
of  measurable  sensitivity  was  calculated  to  be  1143.64  Hz. 
To  verify  this,  the  pitot  static  tube  was  placed  at  the 
center-line  of  the  jet  and  the  Beat  Frequency  Oscillator 

was  traversed  up  to  a frequency  of  3500  Hz.  The  sound 

2 

pressure  level  was  maintained  at  105  dB  re  20  pN/m  . This 
was  measured  with  a BrLiel  and  Kjaer  1/2  inch  (1.27  cm) 
microphone  placed  at  x/D  = 45,  y = 8 i n. (20.32  cm)  and  z = 
8.25  in.  (20.955  cm).  The  microphone  was  connected  to  a 
B & K 2107  Frequency  Analyzer  operating  in  the  linear,  non 
filtering  mode  (Figure  3-9).  By  traversing  the  Beat  Fre- 
quency osc  i 1 1 ator  , the  first  measurable  frequency  at  which 
the  jet  was  most  sensitive  was  found  to  be  1149  Hz.  This 
was  determined  by  watching  the  pressure  fluctuation  in  a 
micro-manometer  as  the  oscillator  settings  were  changed. 
This  value  agreed  with  that  calculated  using  Rockwell's 
(1971)  curve. 
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The  sound  measurements  were  taken  at  three  x/D  sta- 
tions varying  the  lateral  locations.  In  all  the  cases, 
the  sound  pressure  level  was  maintained  at  105  dB,  and  the 
frequency  of  the  acoustic  disturbance  was  held  at  1149  Hz. 
A block  diagram  of  the  equipment  used  to  take  the  sound 
measurements  is  shown  in  Figure  3-10. 


3.  Preliminary  Measurements 
To  determine  how  well  the  jet  agreed  with  published 
results  of  previous  investigators,  a series  of  preliminary 
measurements  were  taken.  Specifically,  the  mean  velocity 
profiles  for  the  undisturbed  and  acoustically  excited  cases 
were  measured  at  three  x/D  stations  and  are  shown  in  Figures 
3-11  and  3-12  respectively.  It  appears  that  Goertler's 
(1942)  equation  is  the  most  applicable  to  the  two  dimen- 
sional jet.  This  equation  takes  the  form 


— = 1 - tanh^  n (3-1  ) 

F 


m 

where  n = (o  y)/x  and  o is  a characteristic  of  the  jet. 
Equation  (3-1)  can  be  written  as: 


2 

tanh  [a  y/x] 


5_  = 1 - tanh^  [-^(b/D)(y/b)] 


(3-2) 
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The  value  for  o can  ho  computed  at  any  x/0  station.  The 
velocity  half  width,  b,  is  defined  as  the  location  where 
the  mean  velocity  is  equal  to  one  li  a 1 f the  c e n t e r - 1 i n e 
velocity  at  that  x/D  station.  Hence,  at  y = b,  = 0.5. 

Solving  equation  (3-2)  for  x/D  = 35  and  y/b  = 1,  the  value 
of  o (for  the  undisturbed  case)  was  found  to  be  9.40, 
and  at  x/U  = 55  and  y/b  -•  1 the  value  of  o was  found  to  be 
9.46.  Hence,  o remained  essentially  constant  at  all  x/D 
stations  in  t ti e similarity  region.  For  t ti e acoustically 
excited  case,  o was  found  to  he  8.28  at  x/D  = 35  and  this 
value  also  remained  constant  at  all  other  x/D  stations 
in  the  si m i 1 a r i t y region. 

Measurements  to  assure  that  the  jet  obeyed  the  prin- 
ciples of  similarity  and  self  preservation  as  discussed  by 
Lin  ( 1 959  ),  Abramovich  ( 1 963  ),  Sc h 1 i c h 1 i nq ( 1 968  ) . Hinze  ( 1 959  ) 

and  Tennekes  and  1.  urn  ley  (1972  ),  were  taken.  One  such 
principle  states  that  the  velocity  half  width  is  directly 
proportional  to  the  distance  downstream  of  the  jet  exit. 

Stated  mathematically,  this  principle  will  take  the  form; 

^ = k,(x/D  - C^)  (3-3) 

where  is  the  widening  rate  of  the  velocity  field  and 
is  the  geometric  virtual  origin  of  the  .iet.  Measurements 
reported  by  investigators  like  Van  der  Hegge  /i.inen  ( 1 957  ), 
Miller  and  Commings  (1957),  Heskestad  (19b5),  Rradbury  (1965), 
Jenkins  ( 1 974  ) and  M w 1 e j ( 1 9 75  ) confirm  the  linear 


41 


relationship  given  in  equation  (3-2).  The  jet  used  by 
the  author  was  the  same  as  that  used  by  Jenkins  (1974  and 
Mulej  (1975).  Since  the  jet  was  built  by  Jenkins  (1974), 
it  was  necessary  for  him  to  take  several  points  to  confirm 
that  the  principles  of  similarity  and  self  preservation 
were  obeyed  by  his  jet.  To  check  if  the  slope  and  inter- 
cept of  this  straight  line  agree  with  those  reported  by 
Jenkins  (1974)  and  Mulej  (1975),  the  author  merely  used 
three  points  to  construct  the  straight  line.  More  points, 
however,  would  offer  less  experimental  scatter.  The  three 
points  used  were  found  to  be  sufficient  (Figure  3-13). 

The  values  for  k^  and  C-j  were  found  to  agree  with  those 
reported.  The  resulting  equations  for  the  jet  are: 


(a)  Undisturbed  Case  (see  Figure  3-13) 


^ = 0.092  (x/D  + 0.7) 


(3.4) 


(b)  Acoustically  Excited  Case  (see  Figure  3-14) 


^ = 0.1  (x/D  + 2.25) 


(3.5) 


The  other  test  for  similarity  involves  the  decay  of 
the  center-line  velocity.  It  states  that  the  square  of 
the  center-line  velocity  is  inversely  proportional  to  the 
distance  downstream  of  the  jet  exit  plane.  This  can  be 
represented  by  the  equation: 
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x/D 


Figure  3-14. 


Jet  Widening  Rate  ; 


Excited. 
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(Ti. 


■D. 


-2 


= (x/0  - C^) 


(3-6) 


where  is  the  decay  rate  of  the  center-line  velocity  and 
is  the  kinematic  virtual  origin  of  the  jet.  Flora  and 
Goldschmidt  (1969)  noted  that  there  was  no  reason  to  expect 
any  relationship  between  the  two  virtual  origins.  The 
undisturbed  and  disturbed  cases  are  shown  in  Figures  3-15 
and  3-16,  respectively.  The  resulting  equations  are; 

(a)  Undisturbed  Case  (Figure  3-15) 


fU 

_n 

lu. 


0.  182  (x/D  - 3.30) 


(3-7) 


(b)  Acoustically  Excited  Case  (Figure  3-16) 


U 


= 0.185  (x/D  - 2.97) 


(3-8) 


4 . C 0 jn  p a r i s 0 n with  o t h e r S h e a^  M o wj^ 

The  mean  velocity  profiles  (see  Figures  3-11  and  3-12) 
agree  with  the  theoretical  results  of  Goertlcr  (1942)  (see 
also  equa t i on ( 3 - 1 ) ) . The  principles  of  similarity  and  self 
preservation  are  satisfied  and  the  empirical  constants  ob- 
tained from  experiments  compare  adequately  with  those  of 
previous  investigators.  The  values  are  tabulated  in  Table 

4 

3-1.  For  the  range  of  Reynolds  number  (0.45  x 10  ^ 

4 ' 

8.0  x 10  ) considered  by  the  investigators,  the  widening 


Table  3-1.  Characteristic  Parameters  of  Plane  Jets  Compared. 
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rate  varied  from  0.085  to  0.110  and  the  decay  rate  of  the 
center-line  velocity  varied  from  0.150  to  0.364.  The  rather 
high  value,  0.364,  was  for  the  case  of  mixing  at  the  jet 
exit.  The  two  dimensionality  of  the  jet  used  by  the  author 
had  been  established  by  Jenkins  (1974)  and  was  not  repeated 
here . 

The  jet  widening  rate  at  a Reynolds  number  of  1.78  x 

4 

10  increased  by  8%  when  an  acoustic  disturbance  was  ap- 
plied. This  percent  increase  was  compared  with  those  ob- 
tained by  Kaiser  (1971).  He  reported  a 21.1%  increase  for 

4 

a Reynolds  number  of  0.45  x 10  and  a 13.7%  increase  for 

4 

a Reynolds  number  of  0.60  x 10  . These  values  confirm 
Rockwell's  (1971)  claim  that  the  Reynolds  number  is  an 
important  parameter  in  determining  the  maximum  sensitivity 
of  jets  to  an  applied  disturbance. 

The  objective  of  taking  these  preliminary  measurements 
was  to  check  if  the  flow  field  is  indeed  a plane  turbulent 
jet.  These  measurements  confirm  the  validity  of  the  jet 
and  hence,  it  can  be  used  as  a valid  flow  field  to  proceed 
with  the  proposed  objectives  of  measuring  inter mittency, 
fold-over  and  crossing  frequency  of  the  turbulent- 
nonturbulent  interface. 
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5 . Detection  of  the  Turbulent  Nonturbulefit  Interface 

a . I ntroducti on 

The  original  objective  of  this  research  was  to  see  if 
a traveling  acoustic  wave  could  change  the  occurrence  of 
fold-over  of  the  interface  of  a plane  turbulent  jet.  It 
was  realized  that  in  order  to  achieve  this  goal  a more 
definite  method  than  that  used  by  Mulej  (1975)  needed  to 
be  developed.  Mulej  (1975)  examined  visually  a continuous 
visicorder  record  of  two  hot  wire  voltages.  As  will  be 
discussed  later,  the  ability  to  measure  fold-over  credibly 
depends  on  how  well  the  tur bu 1 ent -non t u r bu 1 en t interface 
is  detected  for  each  of  the  two  hot  wire  voltages.  Since 
the  hot  wire  probes  used  by  Mulej  (1975)  were  positioned 
at  approximately  0.10  in. (0.254  cm)  apart,  the  difficulty 
of  visually  detecting,  at  all  times,  which  of  the  two 
probes  first  becomes  turbulent  or  nonturbulent  cannot  be 
overemphasized  (see  Figure  3-7).  Hence  the  objective  of 
this  research  expanded  to  include  detecting  the  turbulent- 
nonturbulent  interface  by  a digital  analysis.  This  makes 
it  possible  to  process  a larger  amount  of  data  and  hence  to 
get  a more  meaningful  average. 

Before  presenting  the  details  of  the  actual  method 
used,  a brief  review  of  some  of  the  techniques  used  by 
previous  investigators  will  now  be  given.  The  purpose  of 
this  review  is  to  show  how  arbitrary  the  interface  de- 
tection techniques  are. 
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b.  Review  of  Some  Previous  Techniques 

Several  different  methods  have  been  used  to  separate 
a continuous  signal  originating  from  an  intermittently 
turbulent  flow  into  "turbulent"  and  " non tu rbu 1 en t " zones. 
These  methods  can  be  grouped  into  two.  The  first  method 
is  to  examine  a continuous  record  of  the  signal,  either 
visually,  or  by  computer  analysis  (Kaplan  and  Laufer  (1968), 
LaRue  (1973),  Medley  and  Keffer  (1974),  Spencer  (1970), 
etc.)  when  the  length  of  record  available  for  inspection 
is  considerably  larger  than  the  length  of  a typical  "tur- 
bulent burst."  The  second,  the  analog  approach,  is  char- 
acterized by  a built-in  decision  time  and  a detection 
delay.  This  detection  delay  time  is  set  between  the  actual 
event  and  the  time  at  which  the  decision  is  made.  This 
decision  is  made  on  the  basis  of  the  past  history  of  the 
signal  without  the  benefit  of  knowing  the  future.  Kaplan 
and  Laufer  (1968)  and  Kibens  and  Kovasznay  (1970)  inde- 
pendently came  to  the  conclusion  that  an  important  require- 
ment for  good  detection  is  a knowledge  of  the  future  of  the 
signal.  This  knowledge  can  easily  be  achieved  through  a 
digital  computer  method  since  digitized  signals  can  be 
stored.  The  analog  method  requires  some  time  to  generate 
the  detector  function  (the  i ntermi ttency  function,  I(t)) 
and  consequently,  the  signal  from  other  wires  has  to  be 
delayed  to  account  for  this  time  delay. 

Both  the  digital  and  the  analog  methods  have  been 
proved  credible.  The  choice  between  these  two  methods 
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depends  on  the  cost  of  laboratory  time  versus  computer  time. 
As  observed  by  Kibens  and  Kovasznay  (1970)  and  other  re- 
searchers, a certain  amount  cf  arbitrariness  in  forming  the 
intermittency  function,  I{T),  is  inevitable.  The 
difficulty  in  defining  the  location  of  the  turbulent- 

nonturbulent  interface  has  been  compared  with  that  of  the 
shoreline  of  the  sea  (see  Kibens  and  Kovasznay  (1970)).  The 
location  is  dependent  on  the  mag n 1 f i ca t i o n , or  coarseness  of 
the  method  of  detection.  Depending  on  the  size  of  details 

taken  into  account,  the  shoreline  of  the  sea  could  on  one 
hand  be  considered  relatively  smooth  and  on  the  other  hand, 
full  of  inlets  and  islands. 

The  choice  of  an  input  signal  referred  to  by  Kaplan 
and  Laufer  (1968)  as  the  criterion  function,  is  an  impor- 
tant step  in  the  discrimination  process.  This  signal 
must  offer  the  best  contrast  between  the  turbulent  and  non- 
turbulent  regions.  In  general,  the  criterion  function 
should  be  small  in  the  nonturbulent  regions,  and  appreci- 
ably large  in  the  turbulent  regions. 

Townsend  (1949),  investigating  the  turbulent  wake  of 
a circular  cylinder,  used  the  velocity  time  derivative  as 
the  input  signal.  Corrsin  and  Kistler  (1954)  also  used  the 
velocity  time  derivative  for  their  work  in  the  boundary 
layer,  but  suggested  that  an  improved  detection  would 
result  if  I(t)  were  formed  from  a signal  proportional  to 
the  instantaneous  vorticity  level  in  the  flow.  Fiedler 
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and  Head  (1966)  modified  the  circuit  used  hy  Corrsin  and 
Kistler  (1954)  by  pha se- s h i f t i ng  the  rectified  signal  by  a 
time  which  was  small  compared  to  the  period  between  suc- 
cessive turbulent  bursts,  and  adding  it  to  the  non-shifted 
rec  ti f i ed  signal. 

Kaplan  and  Laufer  (1968)  used  the  short-time  variance 
of  the  time  derivative  of  the  hot  wire  voltage  as  the  cri- 
terion function.  The  variance  was  computed  over  a relatively 
short  time  interval  which  varied  from  2.5  to  8 msec.  The 
computed  variance  was  compared  to  a threshold  level  de- 
rived by  taking  a fraction  of  the  local  maxima  of  the 
criterion  function.  If  the  variance  exceeded  the  threshold 
level,  the  i ntermi ttency  function,  I(t),  was  given  the 
value  of  unity,  otherwise  zero. 

Kibens  and  Kovasznay  (1970)  used  the  time  derivative 
of  a component  of  the  vorticity.  Jenkins  (1974)  also  used 
the  same  input  signal.  LaRue  (1973)  claimed  that  the 
temperature  signal  of  the  wake  of  a heated  rod  offered  a 
good  contrast  between  the  turbulent  and  nonturbulent  zones. 
Antonia  (1972),  in  the  boundary  layer,  used  the  time  de- 
rivative of  the  Reynolds  shear  stress  while  in  the  mixing 
layer.  Wygnanski  and  Fiedler  (1970)  used  the  sum  of  the 
squares  of  the  first  and  second  time  derivatives  of  the 
longitudinal  velocity. 


Most  investigators  agree  that  the  best  criterion 
function  is  that  based  on  the  local  vorticity,  which  is 
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/oro  in  tiu'  turn  t u rim  1 1'li  t regions,  ami  a [i  p r oc  i a t>  1 y large 
i ti  the  Uirtnilont  zones.  However,  the  vorticity  is  rather 
il  it  tic  lilt  to  measure  since,  the  measurement  of  even  one 
component  of  the  vorticity  vector  would  reguire  four 
prt'cisely  positioned,  and  nearly  identical,  tiot-wire  probes, 
because  ol  this  difficulty,  most  ex  per  i iiu' liters  use  other 
physical  variables,  or  a combination  of  them,  as  their 
input  siiinal  evt'ii  though  tiiose  signals  offer  less  contrast 
between  the  turbulent  and  minlurbulent  .'ones.  In  the 
plane  Jet,  for  exaiiifile,  lleskestad  used  tlie  band  passed 

velocity  si)uared,  while  Thomas  (1973)  used  the  sum  of  the 
sguare  of  the  band  (lassed  vi'locity  and  its  derivative  in  the 
plane  turbuli'iit  wake. 

Will'll  an  adeguate  input  signal  is  chosen,  the  next 
im()ortant  steps  involve  the  choice  of  the  detection  cri- 
teria, which  consist  of  two  adjustable  parameters:  the 
threshold  level,  and  the  hold  time,  i . Ihe  threshold 
level,  or  the  gate  level,  can  be  defined  as  that  level  to 
which  the  in[mt  signal,  S(t),  is  compared  to  determine 
whether  the  i n t erm i t t ency  function.  I(t),  would  be  set 
to  unity  or  zero.  In  general,  therefore,  the  i n t e rm i t t ency 
function  is  set  egual  to  unity  when  the  input  signal  is 
egual  to  or  greater  than  the  threshold  level.  Otherwise 
the  i 11 1 erm  i 1 1 ency  function  is  set  ei)ual  to  zero. 

The  other  parameter,  the  hold  time,  i ^ , is  the  length 
of  time  intervals  over  which  the  threshold  criterion  is 
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applied.  This  implies  that  the  i ntermi ttency  function, 
l(t),  is  allowed  to  change  only  when  the  input  signal 
crosses  the  threshold  level  for  a sufficiently  long  time. 

Selection  of  nominal  values  for  tliose  parameters, 
threshold  level  and  hold  time,  is  again  rather  arbitrary. 
Each  investigator  selects  the  nominal  values  best  suited 
to  his  flow  field. 

Kaplan  and  Laufer  (1968)  selected  their  threshold 
level  by  choosing  a fraction  of  the  local  maxima  of  the 
criterion  function  over  a time  interval  they  considered 
reasonable.  The  short  time  variance  which  constituted 
their  criterion  function  had  no  "zero-crossings."  No 
hold  time  was  imposed  on  their  criterion  function  to  de- 
termine the  i ntermi ttency  function.  However,  by  taking 
the  short  time  variance  of  the  time  derivative  of  the 
hot  wire  voltage,  they  were  in  effect  smoothing  the  deriva- 
tive signal.  This  smoothing  was  tantamount  to  applying  a 
hold  time,  , , to  the  square  of  the  derivative  of  the  hot 
wire  voltage. 

Kibens  and  Kovasznay  (1970)  chose  the  nominal  value 
of  the  threshold  level  by  examining  how  the  inter mi ttency 
factor,  'i  , and  the  crossing  frequency,  , varied  with 
different  threshold  levels.  Their  thresliold  was  a 
fixed  constant  determined  by  the  threshold  circuitry, 
but  the  threshold  level  (gain),  which  they  defined  as  the 
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ratio  of  the  rms  of  the  turbulent  part  of  the  input  signal 
S(t),  to  the  fixed  threshold  value,  varied.  They  found 
that  between  certain  threshold  levels,  the  i ntermi ttency 
factor,  and  the  crossing  frequency  approaches  some  asymp- 
totic value.  Within  this  range  of  threshold  levels,  there- 
fore, both  the  i ntermi ttency  factor,  and  the  crossing 
frequency  varied  very  slightly.  Their  final  choice  of  a 
threshold  level  for  a given  location  was  arbitrarily  made 
within  the  range  of  acceptable  threshold  levels. 

Kibens  and  Kovasznay  (1970)  chose  their  hold  time 
differently  from  the  way  they  chose  the  threshold  level. 

The  hold  time  was  related  to  a characteristic  time,  T^j^, 
or  a characteristic  frequency,  These  were  defined 

by: 

5'  = 271 F.  S'  = p-  S'  (3-9) 

ch 

where  5 = (3-10) 

§ is  the  time  derivative  of  the  input  signal,  S ( t ) ; S'  and 
5'  are  the  rms  values  of  S(t)  and  S(t),  respectively.  They 
noted  that  the  choice  of  the  nominal  value  of  the  hold 
time,  T^,  should  be  proportional  to  T^^.  Their  criterion 
for  choosing  was  that  ^ 1.  Specifically,  they 

chose  the  ratio  to  be  4/3.  They  noted  that  in  the  boundary 
layer  in  which  they  were  measuring,  the  characteristic 
frequency  was  330  Hz  at  y/6  = 0.26.  The  characteristic 
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time  was  computed  from  ^ — = 3 msec.  Consequently, 

c h 

they  chose  the  hold  time,  t ^ , to  be  4 msec.  This  value 
was  held  constant  for  the  rest  of  their  measurements. 

LaRue  (1973)  selected  his  threshold  level  in  a manner 
similar  to  Kibens  and  Kovasznay  (1970).  However,  his 
hold  time  was  selected  differently.  He  noted  that  the 
order  of  the  minimum  duration  of  a fully  turbulent  "patch" 
should  be  twice  the  Kolmogoroff  length  transit  time, 
2(si|^/U^),  where  is  the  free  stream  velocity.  He  con- 
tended that  in  order  to  detect  the  smallest  turbulent 
"patches,"  the  hold  time  should  be  chosen  to  be  near  this 
value.  Since  he  analyzed  his  data  digitally,  with  a time 
step  between  digital  samples  of  239  ps,  he  chose  a nominal 
hold  time  of  240  vis,  or  3(s:.|^/U^).  Both  the  threshold 
level  and  the  hold  time  were  subsequently  held  constant  for 
the  rest  of  his  experiment.  This  method  of  choosing  the 
hold  time  was  also  used  by  Hedley  and  Keffer  (1974).  Their 
input  signal,  however, was  a double-component  function  which 
involved  the  derivatives  of  U and  V separately. 

Jenkins  (1974)  followed  a procedure  similar  to  Kibens 
and  Kovasznay  (1970)  in  selecting  his  detection  criteria. 

The  detector  circuitry  and  choice  of  input  signal  were 
similar,  although  the  flow  was  not.  While  Kibens  and  Kovas- 
znay (1970)  were  concerned  with  the  boundary  layers,  Jenkins 
(1974)  considered  a plane  turbulent  jet.  He  selected  a 
threshold  level  by  a method  almost  similar  to  Kbens  and 
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Kovasznay  (1970).  Both  Kibens  and  Kovasznay  (1970)  as 
well  as  LaRue  (1973)  report  that  both  the  crossing  fre- 
quency, F , and  the  hold  time,  lut  should  be  relatively 

I H 

insensitive  to  the  threshold  level  selected  to  be  credible 
in  the  discrimination  process.  Jenkins  (1974)  obtained 
a range  of  threshold  level  over  which  the  crossing  fre- 
quency, , was  relatively  flat.  Within  this  range,  the 
i ntermi t tency  factor,  y,  varied  almost  linearly.  It  seems, 
therefore,  that  he  based  his  choice  of  a threshold  level 
only  on  the  crossing  frequency.  He  noted,  however,  that 
at  his  nominal  threshold  level,  the  i ntermi ttency  factor 
reduced  in  slope  (see  Jenkins  (1974)  Figure  3-7). 

A typical  characteristic  frequency,  F^^  measured  by 
Jenkins  ( 1 974  ) at  a relatively  small  x/D  station  and  close  to  the 
center-line  was  in  the  order  of  615  Hz.  He  defined  the  delay  time,  x^, 
which  was  equal  to  the  hold  time,  by  the  relation; 

Substituting  615  Hz  for  F^^^,  he  obtained  a delay  or  hold 
time  equal  to  1.13  msec.  Sweeping  through  the  remainder 
of  the  flow  field  and  finding  corresponding  characteristic 
frequencies,  he  varied  the  nominal  value  of  the  hold  time 
from  3.60  to  7.5  msec,  depending  on  the  position  in  the 
f 1 ow . 

While  the  investigation  techniques  used  by  the  various 
experimenters  reviewed  were  considerably  different,  their 
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i ntermi ttency  and  crossing  frequency  distribution  func- 
tions as  measured  seemed  to  obey  the  same  expressions. 
Corrsin  and  Kistler  (1954)  established  that  the  probability 
density  function  of  the  interface  position  Y,  was  very 
nearly  Gaussian  and  that  the  distribution  of  Y could  be 
expressed  using  the  error  function.  Specifically,  the 
i ntermi ttency  profile  is  approximated  by  the  expression: 

Y(y)  = i [1  - erf  ] (3-12) 

^ rr  o 

where  Y is  the  mean  position  of  the  interface  and  o is 
its  standard  aeviation.  Y is  defined  as  tiie  location  where 
) = .5.  The  crossing  frequency  distribution  is  approxi- 
mated by  the  expression; 

F^(y)  = exp[-(y  - Y)^/(2o^)]  (3-13) 

where  F^^^  is  the  maximum  crossing  frequency  and  occurs  at 

Y. 


c.  Method  Used 

It  was  pointed  out  in  the  last  section  that  the  original 
objective  of  this  research  was  expanded  to  include  detecting 
the  turbulent- non  turbulent  interface  by  digital  analysis. 

The  review  of  some  of  the  previous  techniques  used  by  other 
experimenters  shows  the  degree  of  a rb i t ra r i nes s in  dis- 
criminating between  turbulent  and  nonturhulcnt  zones.  The 
method  used  here  involves  selecting:  (i)  the  input  signal, 
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(ii)  the  threshold  level,  and  (iii)  the  hold  time.  Each  of 
these  will  be  discussed  fully  here  to  establish  the  credi- 
bility of  the  discrimination  process. 

( i ) The  Input  Signal,  $(t) 

The  square  of  the  time  derivative  of  the  hot  wire 
voltage  was  used  as  the  criterion  function  (input  signal). 

This  is  comparable  to  the  variance  signal  used  by  Kaplan 
and  Laufer  (1968),  the  velocity  time  derivative  used  by 
Townsend  (1949),  Corrsin  and  Kistler  (1955),  and  the  square 
of  the  time  derivative  used  by  Fiedler  and  Head  (1966). 

The  voltages  from  two  hot  wire  probes  placed  0.125  in. 

(0.3175  cm)  apart  were  passed  through  a differentiator 
circuit  (see  Appendix  C).  The  resulting  time  derivatives 
of  these  voltages  were  recorded  on  an  analog  tape  recorder 
at  a speed  of  37.5  in/sec  (95.25  cm/sec).  During  the  play 
back,  at  the  same  speed,  the  signals  were  digitized  simul- 
taneously on  a CDC  1700  Hybrid  computer,  using  the  ap- 
propriate software  (assembly  language).  The  sampling 
rate  was  4000  samples/sec  since  the  highest  frequency  of 
interest  was  2KHz.  Further  explanation  for  this  choice  of 
sampling  rate  is  given  in  Appendix  A.  The  data  were  processed 
on  the  Purdue  University  CDC  6400-6500  dual  MACE  digital 
computer.  10.74  sec  of  data  (42960  samples)  were  analyzed. 
This  was  true  for  all  cases  except  those  in  which  the  hold 
time  and  total  time  of  sampling  were  varied  purposely. 

Details  of  the  digitization  process  are  given  in  Appendix  A. 
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Computer  programs  to  check  the  validity  of  the  digitized 
samples  are  presented  and  explained  in  Appendix  B.  A 
comparison  of  Figures  3-7  and  3-8  justifies  the  use  of  the 
time  derivatives  of  the  hot  wire  voltages  rather  than  the 
direct  voltages  in  attempting  to  discriminate  between  the 
turbulent  and  the  nonturbulent  zones. 


( i i ) The  Threshold  Level 

The  threshold  level  used  in  this  work  is  defined  by 
the  relation: 


C 
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(3-14) 
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where  C is  the  threshold  level  which  varies  with  the  gain, 
2 

P ; S'  is  the  mean  square  of  the  input  signal,  S(t). 

Curves  similar  to  Figure  3-17  were  plotted  to  see  how  the 
i ntermi ttency  and  crossing  frequency  varied  with  different 

it 

threshold  levels  at  the  selected  hold  time,  t^,  of  10  msec. 
As  can  be  seen  from  the  figure,  a range  of  threshold  level 
(p  ^ 0.6)  exists  that  makes  F^  relatively  flat.  It  can 
also  be  seen  that  at  P = 0.6,  the  i n termi ttency  has  changed 
slope  drastically,  and  has  started  to  converge  to  an 
asymptotic  value.  Of  course,  both  curves  will  eventually 
go  to  zero  as  P ® . The  selected  value  of  P in  this 

case  was  0.65.  Hence  the  threshold  level  for  this  par- 

2 2 
ticular  location  was  C = 0.65  S'  . 

The  different  values  of  P used  to  plot  Figure  3-17 

were  obtained  from  the  computer  program.  For  example: 

♦ 

The  hold  time,  is  discussed  in  the  next  section. 

H 


Undisturbed  Jet 
T = 10.74  sec,  = 10  msec 


Selected  Threshold  Level 
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P could  be  initialized  to  0.0  and  a "DO  LOOP"  command  im- 
posed on  the  program.  This  would  cause  P to  take  values 
of  0.2,  0.4  etc.  up  to  say,  1.4,  if  the  increment  chosen 
was  0.2.  For  each  value  of  P,  values  of  y and  could 
be  computed.  Curves  (y  and  F^  versus  P)  similar  to 

Figure  3-17  could  then  be  plotted  to  see  if  either  F^  or  y 

had  exhibited  an  asymptotic  behavior.  By  examining  the 

2 

curves,  P*and  hence  the  threshold  level,  C , could  be 
sel ected . 

The  procedure  explained  above  was  repeated  for  each 
location  in  the  flow  field  examined.  The  selected  value 
of  P and  hence  the  threshold  level  varied  with  location. 

The  variation  of  the  selected  values  of  P with  y/b  can  be 
seen  in  Figures  3-18,  3-19  and  3-20  for  x/0  = 35,  45  and 
55  respectively  and  the  hold  time,  t^,  = 10  msec.  The  same 
procedure  was  also  used  for  the  cases  of  the  acoustically 
excited  jet.  Here  again  the  selected  values  of  P varied 
with  location  as  can  also  be  seen  in  Figures  3-18,  3-19  and 
3-20  for  x/D  = 35  , 45  and  55  and  x = 10  msec.  The  curves 

n 

reveal  that  the  threshold  level  i nc rea ses  nonl i nea r 1 y with 
lateral  location. 

Notice  that  for  each  selected  threshold  level,  values 
of  Y and  F^  were  computed.  In  Figure  3-17,  for  example, 
the  selected  value  of  P at  x/D  = 45,  y/b  = 1.66  is  0.65. 

At  this  value  of  P,  the  corresponding  values  of  y and  F^ 
are  0.446  and  13.04  Hz  respectively.  At  other  locations 
P was  selected  to  give  the  correct  y* 
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different  values  of  y and  were  similarly  obtained. 

Figure  3-21  shows  a curve  of  the  selected  values  of  P versus 
the  corresponding  values  of  y at  x/D  = 35  and  = 10  msec. 

( i i i ) The  Hold  Time 

Whereas  the  threshold  level  was  varied  from  location 
to  location  (as  explained  in  the  last  section),  the  hoi  a 
time  was  held  constant  at  10  msec  for  all  cases.  This 
method  of  varying  the  threshold  level  but  maintaining  a 
fixed  hold  time  is  comparable  to  that  used  by  Kibens  and 
Kovasznay  (1970),  Kaplan  and  Laufer  (1968),  Antonia  (1972), 
and  Hedley  and  Keffer  (1974).  However,  the  actual  hold 
time  used  by  those  investigators  varied.  Fo  example 
Kaplan  and  Laufer  ( 1 968)  used  = 0 sec;  Kibens  and 
Kovasznay  (1970)  used  4 msec  and  Antonia  (1972)  used  0.67 
msec.  LaRue  (1973)  fixed  both  the  threshold  level  and  the 
hold  time.  He  chose  = 240  ys  and  the  threshold  level 
was  constant.  Jenkins  (1974)  chose  a fixed  threshold 

level,  and  varied  the  hold  time  from  3.6  msec  to  7.5  msec. 

LaRue  ( 1 973  ) stated  that  an  analog  ci'-cuit  was  neces- 
sary to  study  the  effect  of  hold  time  on  y and  F^.  He 
noted  that  it  was  not  possible  to  program  into  the  computer 
a hold  time  other  than  the  time  between  digital  samples. 
Based  on  his  argument  a hold  time  of  250  us  was  first 
chosen  in  the  initial  phases  of  this  work.  Figure  3-22 
shows  a typical  curve  (similar  to  Figure  3-17)  obtained 
using  this  value  of  (250  us).  It  appears  that  F^ 
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starts  approaching  an  asymptotic  value  at  P = 1.6.  The 
corresponding  crossing  frequency,  , is  approx ima tely 
equal  to  135  Hz.  This  value  is  considered  to  be  high  for 
a plane  turbulent  jet  based  on  the  work  of  Jenkins  (1974) 
and  Mulej  (1975).  Jenkins  (1974)  reported  a value  of 
about  28  Hz.  Furthermore,  the  corresponding  i ntermi ttency 
factor,  y,  at  this  point  approximately  equals  0.325  (see 
Figure  3-22).  However,  the  value  obtained  by  other  investi- 
gators is  approximately  0.45.  This  meant  that  either  tlie 
hold  time  was  incorrectly  chosen  or  that  somehow  the 
samples  were  not  properly  averaged.  It  was  first  thought 
that  the  total  sampling  time  (10.779  sec)  was  too  short; 
it  was  hence  increased  to  16.1685  sec  (64674  data  points) 
to  check  if  indeed  the  sampling  time  was  too  short. 

Figure  3-23  shows  the  results  for  a longer  sampling  time. 
There  is  no  significant  difference  between  it  and  that 
with  a shorter  sampling  time  (Figure  3-22).  It  became 
obvious  then  that  the  sampling  time  was  indeed  sufficiently 
long  but  that  some  method  had  to  be  employed  to  vary  the 
hold  time  and  hence  select  the  "correct"  value. 

To  that  purpose,  data  were  taken  at  the  center-line 
of  the  jet  where  the  flow  should  be  fully  turbulent.  Using 
this  hold  time  (250  usee),  the  values  of  'i  obtained  were: 

Y = 0.877  for  P = 0.4  and  y = 0.238  for  P = 4.0.  The 


"correct"  hold  time  stiould  give  > = 0.99  for  a reasonable 
range  of  P.  This  range  should  be  chosen  such  that 


Figure  3-23.  y and  F_  vs.  P at  x/D  - 45,  y/b  = 1.66  ; T = 16.168 
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2 

C < V where  V is  the  square  of  the  maximum  voltage  of 
I — mdx  mflx 

the  samples.  This  is  further  explained  later.  A computer 
program  (see  Appendix  B)  to  vary  the  hold  time,  was 

developed.  In  this  program,  the  hold  time  was  varied 

from  4 msec  to  12  msec.  At  each  hold  time,  the  threshold 
2 

level,  C ,was  varied  and  the  corresponding  values  of  y and 
computed.  Curves  of  y versus  P for  the  various  values 
of  hold  time  are  shown  in  Figure  3-24.  From  these  curves, 
the  nominal  hold  time  was  chosen  to  be  10  msec.  This  was 
at  one  x/D  station  (x/D  = 45).  It  was  assumed  that  this 
hold  time  would  apply  at  all  x/D  stations.  To  further 
check  this  value  of  hold  time,  the  procedure  explained 
above  was  repeated  at  x/D  = 45  and  y/b  = 1.61.  At  this 
( lateral  location,  the  value  of  y should  be  approximately 

0.5.  This  value  was  obtained  by  referring  to  the  inter- 
mittency  curve  of  Jenkins  (1974)  who  had  compared  his 
curve  to  those  obtained  by  other  investigators  like  Heske- 
stad  (1965)  and  Bradbury  (1965).  Those  investigators  used 
similar  flow  fields  as  Jenkins  (1974)  and  the  author. 

Curves  for  x^  = 4,  5,  7 and  10  msec  are  shown  in  Figures 
3-25,  3-26,  3-27  and  3-28,  respectively.  It  can  be  seen 
from  these  figures  that  the  only  hold  time,  x^,  that  gives 
the  "correct"  i ntermi  ttency  is  Xl,  = 10  msec.  This  correct 

n 

i ntermi ttency  factor  is  y = 0.5.  This  confirms  that 


1 

\ 
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T u = 10  msec  is  the  "correct"  hold  time.  If  hold  times 
H 

(*) 

less  or  greater  than  10  msec  are  chosen,  "false"  values  of 
the  crossing  frequency  will  be  obtained.  This  second  test 
was  applied  at  two  other  x/D  stations  (x/D  = 35,  55)  using 
four  other  hold  time  values  (t^  = 5,  7.5,  9 and  10  msec). 
Again  it  was  found  that  the  i ntermi ttency  factor  was  equal 
to  0.5  only  at  = 10  msec. 

It  was  observed  from  these  experiments  that  the  neces- 
sary and  sufficient  test  for  selecting  the  "correct"  hold 

time  is  to  examine  what  values  of  > are  obtained  for  all 
2 

values  of  C at  the  center-line  of  the  jet  (y/b  = 0).  It 

2 

is  obvious  that  if  P gets  very  large,  C will  also  get 

2 

very  large  (see  equation  3-14).  If,  for  example,  C gets 

larger  than  the  maximum  value  of  the  samples  being  analyzed, 

the  i ntermi ttency  factor,  y,  will  be  zero  for  all  times. 

The  criterion  for  deciding  what  range  of  values  of  P is 

reasonable  is  discussed  shortly. 

As  explained  in  Appendix  B,  a program  which  compares 

the  magnitude  of  each  sample  and  selects  the  maximum  value, 

2 

V , of  these  samples  is  used.  The  mean  square  S'  is 

in  a X 

2 

also  computed.  The  ratio  of  V , to  S'  is  calculated. 

^ max 

This  ratio  helps  decide  on  the  range  of  reasonable  values 
for  P.  A value  of  P greater  than  20  percent  of  this  is 
possibly  unreasonably  high. 

Figure  3-29  shows  sketches  of  the  steps  taken  to  con- 
struct the  i ntermi ttency  function,  I(t).  The  top  sketch 
Figure  3-29(a)  represents  the  hot  wire  voltage;  3-29{b)  is 


(*)  This  may  bo  duo  to  program  used. 
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the  time  derivative  of  (a).  Both  signals  (a)  and  (b)  are 
obtained  by  analog  method.  The  square  of  (b)  is  shown  in 
(c).  This  is  obtained  by  the  computer  and  is  the  input 
signal  S(t).  By  applying  the  threshold  criterion  alone, 
figure  (d)  can  be  obtained.  This  is  the  intermediate 
signal,  H(t).  By  applying  the  hold  time  criterion  to  (d), 
the  I(t)  function  is  obtained.  Note  that  the  intermediate 
signal,  H(t),  is  actually  not  formed  by  the  computer  as 
programmed.  In  the  program  used,  the  threshold  and  hold 
time  criteria  are  applied  simultaneously  so  that  the  I(t) 
function  is  formed  directly  from  S(t). 

6 . Fold-over 

Fol d-over  (folding)  was  defined  in  Chapter  I as  the 
process  by  which  a region  of  turbulent  fluid  has  non- 
turbulent  fluid  between  it  and  the  fully  turbulent  region, 

A fold-over  region  is  a region  of  turbulent  fluid  which 
has  non-turbulent  fluid  between  it  and  the  fully  turbulent 
region.  This  is  what  Paizis  and  Schwarz  (1974)  called  a 
"fold."  A fold-over  event  occurring  at  the  front  of  the 
interface  is  that  event  whereby  a probe,  further  from  the 
axis  of  the  jet  than  a second  probe  at  the  same  longi- 
tudinal (x-di recti  or)  location,  detects  a turbulent  burst 
before  the  second  probe.  Similarly,  a fold-over  event  oc- 
curring at  the  back  of  the  interface  is  that  event  whereby 
the  probe  further  from  the  jet  axis  remains  in  the  turbu- 
lent region  longer  than  the  probe  closer  to  the  axis  of  the 
jet.  A minimum  of  one  fold-over  event  is  required  to  form 
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a fold-over  region.  However,  a?  many  as  four  fold-over 
events  can  form  a single  fold-over  region.  As  will  be 
discussed  shortly,  back  and  front  fold-over  events  are  also 
called  bac k- f 0 1 d -0 ver  and  front-fold-over,  respectively. 

This  work  did  not  investigate  the  actual  shape  or  width  of 
the  fold-over  region. 

Mulej  (1975)  made  the  first  attempt  to  measure  fold- 
over  events  in  a plane  turbulent  jet.  Earlier  attempts 
were  made  to  quantify  fold-over  in  other  flow  fields.  In 
the  boundary  layer,  for  example,  Paizis  (1972)  and  Paizis 
and  Schwarz  (1974)  observed  that  the  turbulent- 
nonturbulent  interface  exhibited  a significant  amount  of 
fold-over;  (they  used  the  word  "folding"  instead  of  fold- 
over).  They  described  the  interface  position  by  two  single 
valued  functions:  H(x,  z,  t)  and  H*(x,  z,  t).  H(x,  z,  t) 
was  the  total  width  of  turbulent  fluid  at  (x,  z)  at  time 
t.  This  was  measured  by  summing  the  turbulence  detector 
output  (they  used  a linear  array  of  twenty  hot  wire  probes). 
The  second  function,  H*(x,  z,  t),  was  the  width  of  the 
turbulent  fluid  that  was  not  part  of  the  fold-over.  This 
was  obtained  as  the  output  of  the  counting  circuit.  The 
difference  of  those  two  functions  was  a measure  of  the 
total  width  of  the  fold-over  region.  They  measured  the 
mean  value  of  that  quantity  to  be  0.056,  where  6 is  the 
boundary  layer  thickness.  They  noted  that  the  total  w i ' • 
of  the  turbulence  was  1.56  and  that  the  standard  d^v 
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of  the  interface  position  was  approximately  0.36.  Compar- 
ing the  total  width  of  the  fold-over  with  this  standard 
deviation,  they  came  to  the  conclusion  that  this  was  sig- 
nificant. They  reported  that  fold-over  occurred  about  40% 
of  the  time  but  that  fold-over  with  a total  width  equal  to 
the  standard  deviation  occurred  with  a probability  less 
than  0.05. 

LaRue  (1973)  called  fold-over  "overhangs."  His  flow 
field  was  the  wake  of  a heated  rod.  He  also  recognized 
that  due  to  fold-over,  the  turbul ent-nontu rbul ent  inter- 
face was  multivalued.  He  used  two  probes:  one  was  fixed 
3t  y/)i.  = 0.349  where  Z is  a characteristic  length.  The 
position  of  the  other  probe  was  varied.  He  noted  that  the 
number  of  fold-over  events  at  the  downstream  interface  edge 
was  two  to  four  times  that  at  the  upstream  edge.  He  esti- 
mated that  the  fraction  of  time  that  there  was  fold-over 
at  the  downstream  edge  varied  from  0.1  to  0.5%. 

Kohan  (1969)  noticed  in  a wall  jet  that  for  two  trans- 
verse locations  in  the  flow  field,  y^  and  y^  where  y2  > y■^t 
the  flow  was  turbulent  at  y2  and  nonturbulent  at  y.j  less 
than  5%  of  the  time.  Schwarz  (1972)  also  made  that  obser- 
vation in  a wall  jet. 

Preliminary  measurements  of  fold-over  in  a plane  jetwere 
taken  by  Mulej  (1975).  These  measurements  were  made  by  visu- 
ally examining  a continuous  and  simultaneous  visicoder  record 
of  two  hot  wire  voltages.  These  voltages  were  for  two 
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probes  separated  in  the  lateral  direction  by  a small  dis- 
tance. In  effect,  the  turbulent-r.onturbulent  interface 
was  visually  detected  by  him.  His  procedure  to  measure 
fold-over  consisted  in  deciding  which  of  the  two  probes 
first  became  turbulent  or  nonturbu  1 ent . The  hot  wire 
probes  positioned  at  approximately  0.1  in  (0.254  cm)  apart 
displayed  very  similar  voltage  signals  and  hence,  posed  a 
problem  of  visually  deciding  which  probe  first  became 
turbulent  or  nonturbul ent . It  has  been  pointed  out  that 
the  difficulty  of  making  those  judgments  at  all  times 
could  lead  to  an  error  in  the  measurement  technique. 

(J)^  was  defined  as  the  percentage  of  fold-over  in 
the  interface.  Thi  s def i n1 ti on  was  retained  in  this 
work,  but  the  measurement  technique  was  improved.  The 
arrangement  of  the  hot  wire  probes  for  fold-over  measure- 
ments is  shown  in  Figure  3-30.  Figure  3-31  (b)  shows  a 
possible  form  of  the  interface  emphasizing  fold-over.  This 
figure  was  taken  from  Mulej  (1975).  Corresponding  time 
derivatives  of  the  hot  wire  voltages  are  shown  in  Figure 
3-3l(a).  The  main  objective  of  this  research  is  to  measure 
fold-over  by  analyzing  the  digitized  signal  with  a digital 
computer  and  to  see  if  a traveling  acoustic  wave  can  change 
its  occurrence.  No  attempt,  therefore,  was  made  to  model 
the  interface  structure. 

The  computer  program  used  to  study  fold-over  is  pre- 
sented in  Appendix  B.  The  logic  is  explained  here.  From 
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Figure  3-30  It  can  be  seen  that  probe  2 is  closer  to  the 
jet  center-line  than  probe  1.  If  the  structure  of  the 
interface  does  not  have  fold-over,  the  i ntermi ttency  func 
tion,  I ( t ),( proper ly  constructed)  will  always  indicate 
probe  2 turbulent  before  probe  1 and  nonturbulent  after 
probe  1.  This  implies  that  the  width,  ^2%  of  the  inter- 
mittency  function  (random  square  wave)  for  probe  2 will 
always  be  greater  than  that  of  probe  1,  w^  (see  Figure  3-31(a)). 
This  is  because  probe  2 is  closer  to  the  jet  center-line 
where  the  flow  should  be  more  turbulent.  The 
front  and  back  of  the  interface  are  also  shown  in  Figure 
3-31 (b)  and  labeled  A and  B,  respectively,  in  Figure  3-31  (a). 
Three  different  types  of  events  can  occur  in  either  the 
front  or  the  back  of  the  interface.  Two  of  these  types  of 
events  are  defined  as  producing  fold-over  and  the  third 
type  is  defined  as  producing  no-fold-over.  The  types  of 
front  events  can  best  be  explained  with  the  aid  of  a 
sketch.  As  soon  as  the  i ntermi ttency  function  has  been 
formed,  each  sample  point  will  have  a value  of  either  zero 
or  unity.  Recall  that  the  i ntermi ttency  function,  I(t), 
was  defined  by: 


I(t)  = 


1 when  flow  is  turbulent 

0 when  flow  is  nonturbulent 


(3-15) 


The  digitized  signals  from  the  two  hot  wire  probes 
can  be  represented  by  a matrix  of  the  form  [i  x j]  where 
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i * 1,  2 and  j « 1,  2,  3,  42960.  The  number,  42960, 

is  the  total  number  of  samples  processed  per  tape  recorder 
channel.  The  signals  from  the  two  channels  were,  of  course 
digitized  simultaneously.  Figure  3-32  shows  the  three 
types  of  events  that  can  occur  at  the  front  of  the  inter- 
face. Consider  Figure  3-  32(a)  and  the  set  of  points, 

(1,  j)  and  (2,  j).  Both  (1,  j-1)  and  (2,  j-1)  are  zero, 
while  the  points  in  question,  (1,  j)  and  (2,  j),are  both 
equal  to  unity.  This  means  that  both  probes  became  turbu- 
lent simultaneously.  This  type  of  event  is  defined  as 
f rent- f ol d-over  type  1.  Recall  that  if  fold-over  were  not 
possible,  point  (1,  j)  would  have  been  equal  to  zero. 

Recall  also  that  probe  2 is  closer  to  the  jet  center-line 
than  probe  1.  Fold-over  type  2 can  be  explained  similarly. 
This  time,  while  points  (1,  j-1)  and  (2,  j-1)  are  both 
equal  to  zero,  point  (1,  j)  = 1 and  point  (2,  j)  = 0 
(Figure  3-32(b)).  Here  probe  1 becomes  turbulent  before 
probe  two,  thus  producing  front-fold-over  type  2.  The 
third  type  of  event  is  the  expected  event  if  fold-over  were 
not  present.  Here  while  points  (1,  j-1)  and  (2,  j-1)  are 
both  equal  to  zero,  point  (1,  j)  * 0 and  (2,  j)  1 im- 
plying that  probe  2 became  turbulent  before  probe  1 
(Figure  3"32(c)).  Here,  no  fold-over  has  occurred  and  the 
type  of  event  is  defined  as  no-fold-over,  type  5. 

The  back  of  the  interface  also  exhibits  fold-over  be- 
havior. Figure  3-33  shows  the  three  types  of  events  that 
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Figure  3-33.  Back  Events, 


can  occur  at  the  back  of  the  interface.  Consider  Figure 
3-33(a).  Points  (1,  j-1)  and  (2,  j-1)  are  both  equal  to 
1 while  points  (1,  j)  = 0 and  (2,  j)  = 0.  This  means  that 
the  two  probes  become  nonturbulent  at  the  same  time, 
thus  producing  back- f o 1 d-over  type  3.  The  other  type  of 
back-fold-over  occurs  when  both  points  (1,  j-1)  and 
(2,  j-1)  are  equal  to  one  and  points  (1,  j)  and  (2,j)  are 
equal  to  one  and  zero,  respectively  (see  Figure  3-33  (b)). 
This  is  termed  back-fold-over  type  4.  The  final  type  of 
event  occurring  at  the  back  of  the  interface  is  no-back- 
fold-over  type  6.  This  occurs  when  points  (1,  j-1)  and 
(2,  j-1)  are  both  equal  to  one,  and  points  (1,  j)  and 
(2,  j)  are  equal  to  zero  and  one,  respectively  (see  Figure 
3-  33(c)). 

The  computer  program  that  analyzed  fold-over  events 
was  developed  to  make  those  decisions  whenever  a turbulent 
burst  crossed  one  of  the  probes.  There  is  a possibility 
of  a strange  event  occurring.  This  will  be  be  explained 
shortly.  It  is  obvious  that  for  every  front  there  has  to 
be  a back.  If,  for  example,  the  front-fold-over  event 
is  type  1 and  the  corresponding  back-fold-over  is  type  3, 
we  shall  get  a situation  where  the  turbulent  burst  is 
detected  simultaneously  by  the  two  probes  for  both  its  front 
and  back.  This  will  be  a burst  which  at  the  location  of 
crossing  the  probes  has  a perfectly  vertical  front  and 
back.  This  type  of  event  is  called  the  strange  event  and 
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is  considered  no-fold-over  event.  The  computer  program  was 
purposely  written  to  check  for  this  event.  This  strange 
event  was  not  detected  in  all  the  cases  considered,  both 
for  the  disturbed  and  undisturbed  jet.  This  suggests  that 
the  interface  was  either  in  a perfect  fold-over  condi- 
tion or  in  a perfect  no-fold-over  condition.  Fold-over 
types  1 and  3 (i.c.,  the  situations  where  both  probes  be- 
come respectively  turbulent  at  the  same  time,  and  non- 
turbulent  simultaneously)  were  found  to  bo  relatively  few 
in  number.  Fold-over  type  1 was  on  the  order  of  5 percent 
of  type  2.  The  same  was  true  at  the  back. 

In  the  logic  used  to  define  the  six  types  of  events 
occurring  at  the  interface,  it  was  necessary  to  know  what 
the  immediate  past  history  of  the  signal  was.  For  example, 
in  order  to  have  front-fold-over  type  1 at  the  location 
defined  by  the  points  (1,  j)  and  (2,  j),  it  is  necessary, 
but  not  sufficient,  to  have  the  values  of  points  (1,  j-1) 
and  (2,  j-1)  both  equal  to  zero.  The  sufficiency  is  sup- 
plied by  having  points  (1,  j)  and  (2,  j)  both  equal  to 
unity.  Using  mathematical  notation,  front - fol d-over  type  1 
can  be  defined  by  the  relation: 


Type 


(1)  (1.  j-1)  = (2,  j-1)  = 0 
.(2)  (1,  j)  = (2,  j)  = 1 


(3-16) 
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Similarly,  front- fo 1 d-over  type  2 can  be  defined  by: 


Type  2=2  <=> 


'(1)  (1.  j-1)  = (2.  j-1)  = 0 

(2)  (1.  j)  = 1 (3-17) 

(3)  (2.  j)  = 0 


The  rest  of  the  types  of  events  can  be  defined  in  the  same 
way.  (percent  fold-over)  is  defined  by: 

4 

E Type  N Events 

4’i  = X 100  (3-18) 

E Type  M Events 
M=1 


and  fold-over  frequency,  Fg,  is  defined  by: 
4 

E Type  N Events 


where  T is  the  total  time  of  sampling. 

As  was  previously  remarked,  the  accuracy  of  the  fold- 
over  analysis  depends  on  how  well  the  interm ittency  func- 
tion, I(t),  was  formed.  The  logic  used  in  forming  I(t) 
is  presented  in  Appendix  B. 

The  accuracy  in  measuring  fold-over  could  also  be 
affected  by  the  choice  of  Ay  (the  probe  separation).  The 
relationship  of  the  measured  fold-over  percent  to  different 
values  of  Ay  is  shown  in  Figure  3-34.  The  error  associ- 
ated with  the  choice  of  Ay  is  discussed  in  Chapter  V.  It 
should  be  pointed  out  here  that  at  Ay  < 0.0625  in. (0.15875 
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cm)  probe  interference  seems  to  be  present  and  at  Ay  > 
0.1575  in. (0.4  cm)  the  probe  separation  begins  to  exceed 
the  microscale  of  the  flow  and  the  correlation  fails 
(Jenkins  (1974)).  Favre,  Gaviglio  and  Dumas  (1958),  while 
measuring  Space-Time  Correlations  of  velicity  in  a turbu- 
lent boundary  layer,  show  that  no  appreciable  errors  exist 
by  slightly  displacing  one  probe  laterally,  with  respect 
to  the  other.  The  same  observation  was  made  by  Champagne, 
Harris  and  Corrsin  (1970).  From  Figure  3-26,  it  can  be 
seen  that  the  percent  fold-over  is  essentially  constant 
within  a range  of  probe  separation  given  by  0.093  in. 

Ay  _<  0.1  575  in.  (0.236  cm  £ Ay  _<  0.4  cm).  The  final 
value  selected  was  0.125  in. (0.3175  cm).  This  was  chosen 
based  on  Figure  3-34  and  also  for  the  purpose  of  conven- 
ient measurement. 
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CHAPTER  IV.  EXPERIMENTAL  RESULTS 

As  was  already  stated,  the  objectives  of  this  research 
are  three-fold:  (1)  to  detect  the  turbul ent-nonturbul ent 
interface,  (2)  to  measure  fold-over  and  (3)  to  see  the  ef- 
fects of  a plane  traveling  acoustic  wave  on  fold-over, 
crossing  frequency  and  i ntermi ttency . These  objectives 
were  to  be  fulfilled  by  the  use  of  a digital  computer. 

The  results  of  these  measurements  are  presented  in 
this  chapter.  The  i ntermi ttency  and  crossing  frequency 
measurements  were  possible  only  after  a satisfactory  tech- 
nique for  detecting  the  turbul ent-nonturbul ent  interface 
was  developed.  Fold-over  measurements  were  taken  after 
the  i ntermi ttency  function,  I(t),  was  formed.  The  results 
for  the  i ntermi ttency  and  crossing  frequency  with  and 
without  an  acoustic  disturbance  will  be  presented  first. 
Fold-over  results  with  and  without  an  acoustic  disturbance 
will  be  given  next.  These  results  will  be  compared  with 
those  of  other  investigators  in  Chapter  V. 

1 . I ntermi ttency  and  Crossing  Frequency 

All  the  measurements  were  taken  at  a Reynolds  number 
of  1.78  X 10^.  Three  x/D  stations  were  considered  (35, 

45,  55).  The  curve  of  i ntermi ttency  versus  y/b  is  shown 
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in  Figure  4-1.  This  is  the  case  for  the  undisturbed  jet.  A 
few  points  from  Jenkins'  (1974)  curve  are  plotted  to  see  how 
well  the  present  work  compared  with  his  published  work.  It 
can  be  seen  that  the  curve  agrees  with  Jenkins  (1974).  The 
dashed  line  is  an  error  function  curve  fit  of  the  form  given 
in  equation  3-12.  This  equation  can  be  written  in  non- 
dimensional  form: 


T(y')  = i [1  - erf(^  ■ y )]  (4-1) 

/2~  a 

where  y'  = y/b,  Y'  = Y/b  and  o'  = o/b.  The  coefficients 
for  the  curve  are:  Y'  = 1.61,  o'  = 0.408  and  o = 0.670 
in. (1.70  cm).  This  was  calculated  at  x/D  = 35. 

The  i ntermi ttency  profile  for  the  acoustically  dis- 
turbed jet  is  shown  in  Figure  4-2.  Notice  that  while  the 
shape  of  the  curve  remains  unchanged,  the  profile  for  the 
acoustically  excited  jet  has  been  shifted  to  the  left  of 
the  undisturbed  case.  The  coefficients  for  the  curve  fit  at 

x/D  = 35  are:  Y'  = 1.48,  o'  = 0.369,  a = 0.69  in. (1.75  cm). 

The  crossing  frequency  curve  undis- 

turbed case  is  shown  in  Figure  4-3.  This  approximates  a 
Gaussian  curve  of  the  form  given  in  Equation  3-13.  This 
equation  can  also  be  written  in  non-dimensional  form; 


o x/D  = 45 

4 x/D  = 55 
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. exp[.  (4-2) 

Sm  2 o’^ 

where  y'  = y/b,  Y'  = Y/b,  o'  = o/b  and  b is  the  half  width 
of  the  velocity  field  of  the  jet  at  any  x/0  station.  The 
coefficients  for  the  curve  fit  are:  Y'  = 1.66,  o'  = 0.345, 

0 = 0.567  in. (1.44  cm).  A few  points  from  Jenkins'  (1974) 
curve  are  also  shown  in  the  figure.  There  is  a good  agree- 
ment between  the  two.  However  the  value  of  the  maximum 
crossing  frequency,  obtained  by  Jenkins  ( 1 974  ) is  al- 

most double  the  value  obtained  from  the  present  research. 
This  difference  will  be  discussed  further  in  Chapter  V. 

The  crossing  frequency  curve  (F^/F_^^)  for  the 
acoustically  excited  jet  is  shown  in  Figure  4-4.  The  curve 
also  approximates  a Gaussian  curve.  The  coefficients  for 
the  curve  fit  are:  Y'  = 1.43,  o'  = 0.329,  a = 0.613  in. 
(1.56  cm)*.  Notice  that  the  location  Y'  has  decreased.  This 
will  be  discussed  further  in  the  next  chapter.  The  data 
points  used  to  plot  y c>nd  F^  versus  y/b  are  shown  in 
Tables  B-1  and  B-2  for  the  undisturbed  and  acoustically 
disturbed  cases,  respectively. 

In  Figures  4-5,  4-6,  and  4-7,  F^  is  plotted  against 
y/b  for  x/D  = 35,  45,  55,  respectively.  The  acoustically 
disturbed  cases  are  plotted  along  with  the  undisturbed  cases 
to  see  what  effect  the  acoustic  wave  has  on  the  interface 
crossing  frequency.  It  can  be  seen  that  the  crossing 
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frequency  profiles  are  shifted,  and  that  the  maximum  cross- 
ing frequencies  now  occur  at  smaller  values  of  y/b.  Figure 
4-6  shows  that  the  maximum  crossing  frequency  was  slightly 
increased  by  the  acoustic  disturbance.  This  increase  is 
not  noted  in  Figures  4-5  and  4-7  and  this  may  be  due  to 
the  lack  of  experimental  data  points. 

2 . Fold-over  of  the  Interface 

Fold-over  measurements  were  taken  at  the  same  Reynolds 
number  and  x/D  stations  as  those  for  the  i ntermi ttency  and 
crossing  frequency.  The  separation  between  probes  was  0.125 
in.  (0.3175  cm).  A plot  of  the  total  percent  fold-over,  <})^  , 
versus  y/b  for  the  undisturbed  case  is  shown  in  Figure  4-8. 
The  result  shows  that  the  total  percent  fold-over  increases 
with  y/b  and  x/0.  Using  these  figures,  the  range  of  the 
fraction  of  time  during  which  fold-over  events  were  record- 
ed, was  calculated.  Recall  that  42960  points  were  sampled 
and  that  the  total  record  time  was  10.74  sec.  Assuming 
that  42960  fold-over  events  were  recorded,  then  fold-over 
events  would  have  occurred  100%  of  the  time.  It  was  cal- 
culated that  the  fraction  of  time  fold-over  events  were 
recorded  in  this  work  varied  from  0.02%  to  0.42%.  This 
compares  with  0.1  to  0.5%  reported  by  LaRue  (1973)  at  the 
downstream  edge  of  the  interface. 

When  an  acoustic  disturbance  was  applied  to  the  jet, 
the  total  percent  fold-over  was  also  found  to  increase  with 
y/b  and  x/D. 


4-8,  Total  Percent  Fold -Over  vs.  y/b. 
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The  total  percent  fold-over  accounts  for  both  the 
front  and  the  back  fold-over  events.  To  check  whether  the 
fronts  or  the  backs  have  a more  frequent  occurrence  of  fold- 

over,  the  events  are  accordingly  distinguished  (see  Fi cures 

★ * 

4-9a,  4-9b  and  4-9c  ).  These  figures  show  separately  the 
percent  fold-over  for  the  front  and  back  of  the  interface. 
Whereas  Mulej  (1975)  reported  that  fold-over  is  more  likely 
to  be  detected  at  the  front  of  the  interface,  these  figures 
show  that  fold-over  is  actually  more  likely  to  be  detected 
at  the  back  of  the  interface.  Specifically,  Mulej  (1975) 
measured  fold-over  at  11  lateral  locations.  The  percent 
fold-over  for  the  front  was  greater  than  that  for  the  back 
in  10  out  of  the  11  cases  he  considered.  The  present  work 
measured  fold-over  at  15  lateral  locations.  It  can  be  seen 
from  Figures  4-9a,  4-9b,  and  4-9c  that  in  the  average,  the 
percent  fold-over  for  the  back  is  greater  than  that  for  the 
front.  This  is  for  the  undisturbed  case.  Fold-over  mea- 
surements were  taken  at  9 different  lateral  locations  for 
the  acoustically  excited  jet.  It  was  found  that  in  6 out 
of  the  9 cases,  the  percent  fold-over  for  the  back  was 
greater  than  that  for  the  front  (see  Table  B-4).  Further- 
more, examining  each  x/0  station  separately,  it  can  be 
seen  from  Figures  4-9a,  4-9b,  and  4-9c  that  there  is  an 

it 

The  data  show  an  anomalous  behavior  at  y/b  = 1.2.  The 

cause  of  this  is  uncertain. 


no 
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average  percent  increase  for  the  front  over  the  back. 
Specifically,  going  from  the  front  to  the  hack,  there  seems 

to  be  a S%  increase  in  if  ^ at  x/D  = 35  , 55i  increase  at  x/D 

= 45  and  4'i  increase  at  x/D  = 55.  These  results  will  be 
compared  with  those  of  previous  investigators  in  the 
following  chapter. 

To  see  the  effect  of  an  acoustic  disturbance  on  fold- 

over,  Figures  4-10,  4-11,  and  4-12  were  plotted  for  x/D  * 

* 

35,  45,  and  55,  respectively  . It  can  be  seen  that  fold- 
over  occurs  more  frequently  in  the  case  of  acoustic  dis- 
turbance. This  is  more  noticeable  at  x/D  = 35  and  45  than 

at  x/D  = 55.  (The  anomalous  behavior  at  y/b  = 1. 2 is  again 

noticed  in  Figures  4-11  and  4-12). 

It  is  proposed  that  fold- over  could  better  be  repre- 
sented as  a frequency  rather  than  as  a percent.  This  led 
to  a definition  of  fold-over  frequency,  F^,  (see  equation 
3-19).  This  was  defined  as  the  total  fold-over  events 
divided  by  the  total  time  of  sampling.  An  attempt  to 
non-d  i mens  i ona  I i z e F^^,  using  the  maximum  crossing  frequency 

of  the  interface,  F , was  made.  The  plot  of  F,/F 

>m  0 >m 

versus  y/b  is  shown  in  Figure  4-13.  Notice  the  lack  of  a 
universal  trend  for  the  shape  of  the  curve.  The  fold-over 
frequency  was  then  non-d i men s i ona I i z ed  using  the  maximum 

* 

X is  always  measured  from  the  nozzle.  The  shift  of  the 
virtual  origin  is  not  taken  into  account. 


d 
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fold-over  frequency,  The  plot  of  versus  y/b  is 

shown  in  Figure  4-14.  The  results  now  show  a general 
tend.  The  maximum  fold-over  frequencies  were  not  obtained 
very  accurately.  More  data  points  and  larger  samples  are 
needed  to  obtain  better  values.  The  plot  of  F,/F,  versus 
y/b  seems  to  approximate  a Gaussian  curve.  To  this  end, 
a curve  fit  was  attempted  using  equation  4-2.  The  coef- 
ficients for  the  curve  fit  are:  Y'  = 1.58,  o'  * 0.3, 

0 = 0.493  in.  (1.25  cm).  This  is  for  the  undisturbed  jet. 

Figure  4-15  shows  the  plot  of  F^^/F^^  versus  y/b  for 
the  case  of  the  acoustically  excited  jet.  Here  again  the 
curve  is  approximated  by  a Gaussian  curve.  The  coefficients 
for  the  curve  fit  are:  Y'  = 1.48,  o'  = 0.275,  o = 0.513  in. 
(1.302  cm).  To  observe  the  effect  of  the  acoustic  disturb- 
ance on  the  maximum  fold-over  frequencies.  Figures  4-16, 

4-17  and  4-18  were  plotted  for  x/0  = 35,  35,  and  55, 
respectively.  The  results  show  the  possibility  of  a slight 
increase  in  the  maximum  fold-over  frequency.  It  is  clear 
that  more  data  points  are  required  to  study  the  actual 
shape  of  the  curves.  However,  it  was  not  the  objective  of 
this  research  to  do  this. 

In  the  intermitten cy,  crossing  frequency  and  fold-over 
frequency  curves  just  reported,  the  mean  interface  location 
Y_^  was  not  identical  with  the  location,  Y^ , where  the  maxi- 
mum crossing  frequency,  occurred.  Recall  that  the 

mean  interface  location  was  defined  arbitrarily  as  the 
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Figure  4-17.  Fold-Over  Frequency  vs.  y/b  at  x/0  = 45. 
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location  where  > = 0.5,  Ideally  the  two  locations  (Y  and 
) would  be  identical.  However,  experimenters  have  re- 
ported slight  variations  in  those  values  (LaRue,  1973; 
Kibens  and  Kovasznay,  1969,  1970;  Jenkins,  1974,  and  Mulej, 
1975).  LaRue  (1973)  for  example  obtained  Y^  = 4.89  cm, 

Y^  = 4.68.  This  is  a variation  of  about  4.3%.  In  the 
present  work,  a 3%  variation  was  obtained.  This  is  well 
within  the  experimental  scatter  and  not  much  should  be 
concluded  from  it. 

The  cause  of  the  increase  in  the  maximum  interface 
crossing  frequency  observed  in  Figure  4-17  is  not  certain. 
It  may  be  due  to  the  positioning  of  the  exponential  horn. 
Recall  that  the  axis  of  the  horn  was  fixed  at  x/D  = 45  for 
all  the  acoustically  excited  cases.  Note  that  the  curves 
shown  in  Figures  4-16,  4-17,  and  4-18  were  not  approxi- 
mated by  a Gaussian  curve  (the  nond  i men s i ona 1 frequencies 
were ) . 

The  anomalous  behavior  at  y/b  = 1,2  is  again  noticed 
in  Figure  4-18.  The  choice  of  in  this  figure  was  made 
to  agree  with  those  obtained  from  Figures  4-16  and  4-17, 
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CHAPTER  V.  DISCUSSION  OF  RESULTS 

This  work  has  been  purely  experimental.  It  has  ex- 
plored digital  methods  to  analyze  the  data.  Although 
analog  methods  are  equally  acceptable,  the  digital  computer 
methods  have  the  advantage  of  using  less  instrumentation 
that  would  otherwise  contribute  to  error.  Digital  methods 
also  permit  repeated  analysis  with  the  same  base  data. 

This  discussion  will  be  divided  into  four  parts:  (1) 
the  results  of  the  measurements  presented  in  the  previous 
chapter  will  be  compared  with  those  of  previous  experi- 
menters. (2)  Possible  sources  and  estimates  of  error  will 
be  presented  following  comparison  with  previous  work.  (3) 

A discussion  on  entrainment  and  the  interface  slope  will 
follow.  (4)  The  effects  of  an  acoustic  distrubance  on  a 
plane  turbulent  jet  will  be  discussed. 

1 . Compa r i son  with  Previous  Work 
a.  The  Undisturbed  Jet 

The  i ntermi ttency  curve  was  shown  in  Figure  4-1.  The 
results  of  this  curve  were  compared  with  Jenkins  (1974)  and 
were  found  to  agree.  Jenkins  (1974)  in  turn  had  compared 
his  results  with  those  obtained  by  Heskestad  (1965)  and 
Bradbury  (1965)  and  found  agreement.  It  was  pointed  out 
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in  Chapter  jythat  the  shape  of  this  curve  is  approximated 
by  an  error  function  of  the  form  given  in  egua  t i on  (4  - 1 ). 

The  coefficients  of  the  curve  fit  were  also  given.  Libby 
(1975)  proposed  a model  equation  which  described  the 
intermittent  turbulent  flow.  Assuming  that  the  uncondi- 
tioned flow  was  known,  he  showed  that  the  properties  of 
the  conditioned  mean  velocity  and  of  the  i nt erm i t t ency 
could  be  predicted.  He  compared  predicted  and  experimental 
quantities  for  two  mixing  layers  and  a simulated  boundary 
layer  and  found  the  agreement  sa  t i s fa-' tory  . This  model 
could  also  be  applied  to  a plane  jet. 

The  crossing  frequency  was  compared  with  Oonkins  (1974) 
and  the  plot  of  F /F  vs.  y/b  agreed.  However,  as  was 
previously  noted,  the  values  of  the  maximum  crossing  fre- 
quency, obtained  in  this  research  did  not  compare  fav- 

orably w i t li  Jenkins'  ( 1 974  ),  even  though  F / F vs.  y/b 
did  agree^*l  No  other  paper  has  been  published  that  gives 
values  for  the  maximum  crossing  frequency  of  the  interface 
for  a two  dimensional  plane  turbulent  jot.  The  values 
obtained  in  this  research  are  shown  compared  with  Jenkins' 

( 1 9 74  ) in  Table  5 - 1 . Similar  values  o ti  t a i n o d for  other 
shear  flows  are  also  shown  in  Table  5-1.  The  non- 
dimensional  variance  obtained  from  an  error  function  curve 

fit  is  denoted  by  e /Y  . The  other  two  uoiul  i mens  i ona  1 

Y Y 

variances,  crossing  and  fold- 

over  frequencies,  respectively. 


(*)  This  could  be  due  to  the  different  hold  times  used. 


Table  5-1.  Maximum  Interface  Crossing  Frequency  and  Parameter  for  the  Curve  Fit. 
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The  values  of  F „ obtained  by  Jenkins  (1974)  are 
approximately  equal  to  that  obtained  by  Kibens  and 
Kovasznay  (1970)  in  the  boundary  layer.  There  is  no  rea- 
son to  expect  the  values  to  agree.  In  fact,  since  the 
flows  are  entirely  different,  it  seems  reasonable  to 
expect  a difference  in  the  maximum  interface  crossing  fre- 
quency. Besides,  the  i ntermi ttency  curves  do  not  match 
even  though  they  obey  the  same  error  function.  The  constants 
in  equations  (3-12)  and  (3-13)  are  different  for  different 
flow  fields  since  o and  Y are  different.  Hence,  it  would 
not  be  surprising  if  the  values  of  the  maximum  crossing  fre- 
quency F , would  be  indeed  different. 

The  values  of  the  maximum  crossing  frequency  obtained 
by  Jenkins  (1974)  were  used  by  Mulej  (1975)  in  his  final 
presentation  of  the  determined  lateral  velocity  of  the 
interface.  He  shows  that  the  lateral  velocity  of  the 
interface  should  take  the  form: 


(5-1) 


Equation  (5-1)is  Mulej's  (1975)  equation  6.  This  equation 
can  be  rewritten  to  give: 


From  equation  (5-2)  it  can  be  seen  that  the  sign  of  de- 
pends ci.  3Y/3y.  Mulej  ( 1975  ) plotted  b(9Y/9y)  vs.  y/b. 
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This  is  actually  a plot  of  the  magnitude  of  b(3Y/dy)  vs. 
y/b.  From  the  equation  for  y (equation  3-12)  and  its 
shape  shown  in  Figure  4-1,  it  is  obvious  that  the  slope  of 
y,  given  by  either  zero  or  negative  all  the 

time;  i.e.,  3Y/3(y/b)  < 0.  Applying  this  condition  to 
equation  5-2  it  can  be  seen  that  will  always  be  posi- 
tive. 

Mulej  (1975)  measured  the  lateral  velocity  of  the 
interface  in  three  ways.  He  defined  as  a measure  of 

the  interface  velocity  neglecting  the  effects  of  fold-over 
and  V^2  velocity  of  the  interface  considering  the 

effects  of  fold-over.  He  claims  that  V^2  takes  no  account 
of  the  direction  of  the  velocity  of  the  interface  and  thus 
was  not  a good  representation  of  the  velocity.  He  finally 
defined  a third  velocity  of  the  interface,  as  that 

velocity  which  takes  into  account  the  direction  of  the 
velocity  of  the  interface.  It  has  been  argued  earlier  in 
this  chapter  that  the  lateral  velocity  of  the  interface, 
as  defined  by  the  equation  (5-2)  cannot  be  negative.  Hence 
Vi3  does  not  seem  to  be  the  best  measure  of  the  velocity  of 
the  interface.  It  was,  however,  not  the  objective  of  this 
research  to  measure  the  velocity  of  the  interface.  The 
measurements  of  by  Mulej  (1975)  offer  a means  of  com- 
paring the  maximum  interface  crossing  frequency,  ob- 

tained in  this  research  to  that  obtained  by  Jenkins  (1974). 

Since  V.  was  defined  as  a function  of  F it  is  possible  to 

] > 

This  applies  to  the  average  value  of  V^.  The  instantaneous 
value  might  be  negative. 
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compare  the  measured  with  the  predicted  value  using 

Jenkins'  (1974)  and  that  obtained  from  the  present  work. 

Since  F = F at  v = 0.5  the  comparison  has  been  made  at 
Y Y f** 

Y = 0.5  at  three  x/D  stations.  The  results  are  tabulated  in 
Table  5-2.  Notice  that  V^./b  computed  with  F^^^,  obtained 
from  the  present  research,  agrees  with  the  measured  values 
V^^/b  and  much  better  than  that  computed  using 

Jenkins'  ( 1 974  ) values.  However  agrees  better  with 

V^/b  computed  using  Jenkins'  (1974)  F^^, 


Since  takes 


no  account  of  fold-over,  and  defines  the  velocity  of 

the  interface  as  both  positive  and  negative,  it  is  assumed 

here  that  both  defintions  could  give  incorrect  values 

of  the  velocity  of  the  interface.  It  is  suggested  here 

that  a reasonable  definition  of  the  lateral  velocity 

of  the  interface,  as  measured  by  Muiej  (1975),  is  that 

given  by  V^2'  relative  good  agreement  between  the 

measured  velocity,  V^2*  predicted,  V^,  using  the 

present  F suggests  that  the  values  of  F obtained  from 
Y't’  Y'’’ 

this  research  are  probably  the  "correct"  values. 

Both  Kibens  and  Kovasznay  (1970)  and  LaRue  (1973) 
report  that  F^  is  very  sensitive  to  the  variations  in  the 
technique  used  to  construct  the  i n term i ttency  function, 

I(t).  To  show  how  F varies  with  the  hold  time,  Tu.  this 
work.  Figure  5-1  is  shown.  Notice  that  — 5 msec,  F 


31  sec“\ 
Chapter  III 


Ym 


The  method  used  to  choose  t was  presented  in 

. Tm  = 5 msec  was  found  to  be  unacceptable. 

H 


Table  5-2.  V.  Obtained  Using  F Obtained  by  Jenkins  (1974)  and  Present  Work. 
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Undisturbed  Jet 


J 


( i } Fol d-Over 


The  only  reported  measure  of  fold-over  in  a plane 
turbulent  jet  is  that  given  by  Mulej  (1975).  Other  mea- 
surements of  fold-over  in  other  shear  flows  were  discussed 
in  Chapter  III.  It  was  pointed  out  earlier  that  the 
method  used  by  Mulej  (1975)  could  lead  to  erroneous  mea- 
surements of  fold-over.  It  has  been  stated  that  fold-over 
can  be  measured  only  after  the  i ntermi ttency  function, 

I(t),  has  been  formed.  Whereas,  Mulej  (1975)  reported  that 
fold-over  decreased  with  y/b,  this  research  shows  that  it 
increases  with  y/b  instead.  There  is  no  theory  to  sup- 
port this  claim.  However,  it  is  shown  here  that  a de- 
creasing function  can  be  obtained  by  using  incorrect 
threshold  levels.  These  values  uf  the  threshold  level 
would  give  incorrect  values  of  y and  F^.  Figure  5-2  shows 
the  two  trends.  The  data  points  used  to  plot  those  curves 
are  shown  in  Table  5-3.  Notice  that  different  values  of 
(percent  fold-over  accounting  for  both  the  fronts  and 
the  backs  of  the  interface)  can  be  obtained  by  different 
settings  of  P*  Only  three  locations  were  checked  and  these 
were  found  to  be  sufficient  to  show  the  trend  (see  Figure 
5-2)  . 

The  fold-over  frequency,  F^  , reported  in  this  re- 
search appears  to  be  unique  in  the  sense  that  it  has  not 
been  found  to  have  been  reported  by  any  other  investigator. 
This  frequency  seems  to  be  a good  measure  of  fold-over. 

*These  values  of  P were  selected  purposely  to  exhibit  the 
"wrong"  trend. 
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Figure  5-2.  Percent  Fold-Over  vs.  y/b. 
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Table  5-3.  Comparing  Mulej's  (1975)  'P.:  - 10  msec. 
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The  plot  of  F-/F  vs.  y/b  appears  to  approximate  a Gaussian  ^ 

I U V ill 

curve  of  the  form  given  in  equation  4-2. 

It  was  reported  earlier  that  Paizis  and  Schwarz  (1974) 
measured  the  total  width  of  the  fold-over.  Fold-over  widths 
were  not  measured  in  this  research.  LaRue  (1973),  and  also 
LaRue  and  Libby  (1974)  measured  the  fraction  of  time  that 
there  is  fold-over  at  the  downstream  edge  of  the  interface 
of  the  wake  of  a heated  rod.  They  also  reported  that  the 
number  of  fold-over  events  at  the  downstream  interface  edge 
was  two  to  four  times  the  value  measured  at  the  upstream 
interface  edge.  This  means  that  fold-over  is  more  likely 
to  be  detected  at  the  downstream  edge  of  the  interface 
than  at  the  upstream.  The  downstream  interface  edge  of  a 
( wake  corresponds  to  the  back  of  the  interface  in  a jet  and 

the  upstream  corresponds  to  the  front.  This  work  also 
reports  that  fold-over  is  more  frequent  at  the  back  than  at 
the  front  of  the  interface,  but  not  as  much  as  reported  by 

LaRue  and  Libby  (1974), 

Kohan  (1969)  and  Paizis  (1974)  reported  that  in  the 
plane  wall  jet  in  which  they  were  measuring,  fold-over 
events  occurred  less  than  55.  of  the  time.  LaRue  (1973) 
and  LaRue  and  Libby  (1974)  reported  that  fold-over  events 
occurred  at  the  downstream  edge  of  the  interface  between 
0.1  to  0.5«  of  the  time.  In  the  present  work,  fold-over 
events  were  found  to  occur  between  0.2  to  0.425.  of  the  time. 

This  seems  to  show  agreement  with  previous  workers. 


( 
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b.  The  Acoustically  Excited  Jet 
It  was  mentioned  in  Chapter  II  that  no  previous  work 
has  been  reported  concerning  the  effects  of  an  acoustic 
disturbance  on  i n t erm i t tency  , fold-over  and  crossing  fre- 
quency of  the  tu r bu 1 en t -nont u r bu 1 en t interface  of  a shear 
flow.  These  measurements  appear  to  be  unique  and  hence 
are  not  compared  with  other  results. 

2 . Possible  Sources  of  Error  and  Estimates 

a.  Threshold  Level  and  Hold  Time 
The  arbitrariness  used  in  selecting  the  criteria  for 
detecting  the  turbulent-nonturbulent  interface  could  lead 
to  some  error  in  constructing  the  i ntermi ttency  function, 
I(t).  The  method  used  to  select  the  hold  time  was  dis- 
cussed in  Chapter  III.  It  was  suggested  earlier  that  the 
error  involved  in  selecting  was  small.  This  error  is 
estimated  at  ^ 0.5*  of  1^.  However  the  threshold  level 
presented  some  problems.  A range  of  threshold  level 

existed  that  gave  relatively  constant  F values.  However, 

Y 

the  problem  was  in  deciding  what  value  of  threshold  level 
to  choose  within  the  acceptable  range.  The  calibration 
curve.  Figure  3-28,  showed  that  the  selected  threshold  level 
was  close  to  the  start  of  the  acceptable  range.  Within 
this  acceptable  range,  y varied  by  about  21%  and  F^  by  3%. 
However,  when  the  threshold  level  was  chosen  close  to  the 
start  of  the  acceptable  range  of  F^ , the  i ntermi ttency , y, 
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varied  by  3'i  and  by  0.5".  This  method  of  choosing  the 
threshold  level  dose  to  the  start  of  the  acceptable  range 
of  was  adopted  for  the  rest  of  the  experiment.  The  un- 
certainty in  the  present  dat<\  (based  on  the  estimated  error 
in  constructing  l{t))  is  believed  to  bo  -3%. 

b.  Probe  Separation,  Ay 

The  measurement  of  Ay  is  accurate  to  about  0.01  in. 
(0.02 5 A cm).  It  was  discussed  in  Chapter  111  liow  this 
separation  was  achieved  and  how  Ay  was  chosen  for  fold- 
over  measurements.  The  curve  of  vs.  Ay  was  shown  in 
Figure  3-34.  Using  the  value  of  Ay  (0.125  in.  or  0.3175 
cm),  the  error  in  measuring  percent  fold-over  was  found 
to  be  -4*.  However,  bettor  results  might  have  been  ob- 
tained by  takii\g  measurements  with  different  Ay  and 
ex  trapo  I a t i tig  the  results  to  zero  separation. 

c . Errors  Due  to  0 i g i t a 1 Analysis 

These  sources  of  error  are  discussed  in  Appendix  A. 
Some  of  t h 0 s e e r r o r s c o u 1 d be  due  to  ( 1 ) t ti e s a m p 1 i n g rate, 
(2)  total  time  of  sampling,  (3)  aliasing.  Other  errors 
could  be  due  to  the  calibration  of  the  tape  recorder  and 
the  clioice  of  tiie  differentiator  circuit  parameters.  The 
differentiator  circuit  is  discussed  in  Appendix  C. 

d . Cost  E s t i ma  te 

A cost  estimate  of  using  the  digital  computer  is  made 
hero.  In  Chapter  III,  the  method  of  selecting  the  thresh- 
old level  was  discussed.  Recall  that  in  ordi'i’  to  select 
t ti e "correct"  P , different  v .i  1 u e s of  P were  u s (' d a n d curves 


of  Y and  F versus  P were  drawn.  Up  to  33  different  values 
Y 

of  P had  been  used  to  plot  one  of  those  curves  (however, 

fewer  points  wore  used  in  many  cases).  This  means  that  the 

computer  program,  Thold  10  (see  Appendix  B)  iterates  33 

times.  Of  course,  this  was  not  done  at  one  time.  Ten 

values  of  P could  be  selected  at  first;  then  y -Jiid  F could 

Y 

be  plotted  against  P.  The  shape  of  the  F^  curve  would 
decide  whether  to  select  more  values  of  P and  where  to 
select  them.  The  cost  for  each  iteration  is  approximately 
$2.00.  Thirty-three  iterations  would  cost  $66.00.  This 
is  for  one  y/b  location.  Fifteen  locations  were  examined 
for  the  fold-over  measurements  for  the  undisturbed  jet; 
nine  locations  were  considered  for  the  acoustically  excited 
jet.  This  gives  a total  of  24  locations  or  approximately 
$1584.00.  After  the  "correct"  P is  chosen,  the  fold-over 
program,  FOLD,  is  used  to  compute  the  fold-over  events. 

Each  fold-over  program  cost  approximately  $2.20.  This  gives 
a total  of  about  $52.80.  The  cost  of  using  the  other  two 
programs,  SKIP3  and  MSQ,  is  approximately  $200.00.  The 
total  computer  cost  is  approximated  at  $1800.00.  This 
does  not  include  the  cost  encountered  when  each  program  was 
being  developed. 

3 . E n t r a i n m e n t_  and  Interface  Slope 
The  phenomenon  of  entrainment  is  a feature  of  all  free 
turbu’ent  shear  flows.  The  entrainment  process  is  that 
whereby  ambient  nonturbulent  fluid  is  being  assimilated 


by  the  turbulence.  The  exact  process  of  this  mechanism  is 
still  not  fully  understood.  However,  it  is  generally  accept- 
ed that  the  originally  irrotational  fluid  acquires  vorticity 
by  viscous  diffusion  and  this  vorticity  is  amplified  by  the 
rate-of-strain  field.  Paizis  and  Schwarz  (1975)  argue  that 
this  process  must  lead  to  a sharp  interface  which  separates 
the  turbulent  and  the  nonturbulent  regions.  They  defined 
the  rate  of  entrainment  of  ambient  fluid  across  a turbulent 
interface  as  the  mean  rate  of  increase  of  turbulent  fluid  in 
the  flow  direction.  The  flow  field  considered  in  their  ex- 
periment was  a two  dimensional  wall  jet.  They  assumed 
that  the  interface  is  a continuous  surface  whose  location 
is  a random  function  of  space  and  time. 

Many  researchers  (Kaplan  and  Laufer,  1967;  Kohan, 

1969;  Townsend,  1970;  LaRue,  1973;  Davies  and  Yule,  1974; 
Paizis  and  Schwarz,  1974;  Jenkins,  1974;  Mulej,  1975;. and 
others)  agree  that  fold-over  of  the  interface  is  an  impor- 
tant part  of  the  entrainment  mechanism.  Furthermore, 

LaRue  (1973)  argues  that  fold-over  events  are  intimately 
related  to  the  mean  interface  slope.  His  measurements  show 
that  the  mean  slope  at  the  downstream  (back)  interface 
edge  is  steeper  than  the  mean  slope  at  the  upstream  (front) 
interface  edge.  He  attributed  this  to  the  larger  number 
of  fold-over  events  that  were  observed  at  the  back  of  the 


interface.  Mean  interface  slope  measurements  were  not 
taken  in  this  work.  However,  from  the  measurements  taken. 
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j the  sign  of  the  interface  can  be  inferred.  Whenever  a 

fold-over  event  (fold-over  type  2)  occurs  at  the  front  of 
the  interface,  the  interface  slope  is  negative  at  that 
point.  On  the  other  hand  the  slope  of  the  interface  is 
positive  when  a fold-over  event  (fold-over  type  4)  occurs 
at  the  back  of  the  interface.  Hence  the  total  number  of 
fold-over  events  type  2 is  a measure  of  the  total  number 
of  times  the  front  of  the  interface  is  negative  at  that 
location.  Similarly,  the  total  number  of  fold-over  events 
type  4 is  a measure  of  the  total  number  of  times  the  back 
of  the  interface  is  positive  at  the  location  considered. 
Finally,  the  total  number  of  fold-over  events  type  1 and 
3 is  a measure  of  total  number  of  times  the  slope  of  the 
( interface  is  infinite  at  that  location.  Recall  that  the 

types  of  fold-over  events  were  defined  in  Chapter  III. 

4 . Effects  of  an  Acoustic  Disturbance 
The  application  of  an  acoustic  disturbance  at  a fre- 
quency of  measurable  sensitivity  seems  to  have  two  effects 
on  the  flow:  (1)  the  location  at  which  > = 0.5  (T)  is  re- 
duced for  the  i ntermi ttency , crossing  frequency  and  fold- 
over  frequency  curves.  Recall  that  the  widening  rate  of 

the  jet  increased  by  8.7«.  It  appears  that  the  reduction  in 
the  mean  interface  locations  is  of  the  same  order  as  the 
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increase  in  the  widening  rate  of  the  jet.  (2)  The  other  ef- 
fect of  an  acoustic  disturbance  on  the  flow  in  the  increase 
of  the  maximum  crossing  frequency  and  fold-over  frequency. 
Although  this  increase  could  not  be  inferred  at  x/D  = 35  and 
55,  it  was  clearly  noticed  that  x/D  = 45.  Only  three  loca- 
tions (and  hence  three  data  points)  were  examined  for  each 
x/D  station.  It  is  obvious  that  more  data  points  will  be 
needed  to  establish  this  effect  at  all  x/D  stations.  At  x/D 
= 45  where  the  increased  was  noticed,  the  maximum  crossing 
frequency,  increase  by  5.3*  and  the  fold-over  frequency 

increased  by  47.8%.  The  total  percent  fold-over  was  also 
increased  by  the  acoustic  disturbance.  These  increases  are 
approximately:  9%,  4%,  and  6%  at  x/D  = 35,  45,  and  55%, 
respectively.  These  increases  are  of  the  order  of  the  in- 
crease in  the  widening  rate  of  the  jet. 

The  mechanism  by  which  an  acoustic  wave  interacts  with 
a plane  turbulent  jet  is  still  unknown.  Simcox  (1969)  pro- 
posed that  the  mechanism  was  a direct  interaction  of  the 
acoustic  wave  with  the  turbulence.  Hussain  and  Zaman  (1975) 
seem  to  believe  that  the  interaction  is  at  the  initial 
region  of  the  jet.  (Recall  that  their  acoustic  excitation 
was  at  the  plenum  chamber.)  The  present  work  was  conducted 
entirely  in  the  fully  developed  region  (and  in  the  similar- 
ity region).  The  location  and/or  mechanism  of  the  inter- 
action cannot  be  inferred  from  this  work.  Besides,  since 
the  period  between  turbulent  bursts  were  not  measured,  it 
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is  difficult  to  infer  whether  the  period  between  bursts 
was  more  structured  in  the  presence  of  an  acoustic  excita- 
tion than  in  its  absence.  A study  of  this  interaction  is 
presently  being  conducted  by  Chambers  (1976).  The  present 
research  is  related  to  his  work. 
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CHAPTER  VI.  CONCLUSIONS  AND  SUGGESTIONS 
FOR  FUTURE  STUDY 

The  following  conclusions  are  made  based  on  this  work: 
The  FORTRAN  IV  program  developed  for  this  research  sat- 
isfactorily detected  the  t u r bu  1 en t- non t ur bu 1 en t inter- 
face. 

The  sign  of  the  slope  of  the  interface  can  be  deduced 
from  the  types  of  fold-over  events. 

The  maximum  fold-over  frequencies  F^^  were  found  to  be 
less  than  the  maximum  interface  crossing  frequencies, 

^ym  ■ 

The  presence  of  fold-over,  especially  at  the  back  of 
the  interface  edge,  is  consistent  with  the  hypothesis 
that  one  of  the  entrainment  mechanisms  is  large  scale 
"engulfment"  of  the  nonturbulent  ambient  fluid  by  the 
turbulent  fluid.  However,  the  presence  of  fold-over  is 
inconsistent  with  the  Gaussian  interface  model.  This 
model  assumes  that  there  are  no  fold-over  events  at  the 
fronts  and  backs  of  the  interface  (LaRue  (1973)). 

The  total  percent  fold-over  increased  with  the  lateral 
location  (y/b).  Both  the  total  percent  fold-over  and 
the  number  of  fold-over  events  increased  as  the  longi- 
tudinal distance  (x/D)  increased.  This  is  consistent 
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^ with  the  suggestion  (Phillips,  1972),  that  the  slower 

moving  eddies  are  responsible  fo*'  the  entrainment  of 
ambient  fluid.  A larger  number  of  these  eddies  are 
more  likely  to  be  found  further  from  the  mouth  of  the 
jet  than  close  to  it. 

6.  Fold-over  events  were  more  likely  to  be  detected  at  the 
backs  rather  than  the  fronts  of  the  interface. 

7.  An  acoustic  disturbance  at  a frequency  corresponding  to 
the  first  sensitive  frequency  was  found  to  reduce  the 
mean  interface  locations,  for  the  i ntermi ttency , crossing 
frequency  and  fold-over  frequency  curves.  The  percent 
reductions  were  on  the  order  of  the  percent  increase  in 
the  widening  rate  of  the  jet.  The  total  percent  fold- 

^ over  also  increased  and  was  also  on  the  order  of  the 

percent  increase  in  the  widening  rate  of  the  jet.  The 

I 

; maximum  values  of  the  crossing  frequency  and  fold-over 

j 

I frequency  were  also  found  to  increase  (observed  only  at 

3 x/D  = 45  where  the  axis  of  the  horn  was  positioned). 

8.  In  detecting  the  tu r bu 1 en t -non tu r bu 1 en t interface  by 
this  method,  it  was  convenient  to  select  a fixed  hold 
time  and  vary  the  threshold  level.  However,  the  deci- 

Ision  of  the  detection  criteria  is  left  to  the 

researcher. 
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The  following  are  suggestions  for  future  study. 

1.  The  times  between  turbulent  bursts  were  not  computed  in 
the  work  just  completed.  The  computer  program  can  be 
adjusted  easily  to  measure  this  time.  The  number  of 
points  between  the  (1-0)  and  the  following  (0-1)  level 
changes  can  be  computed.  Dividing  this  number  by  the 
sampling  rate,  the  required  time  between  bursts  can  be 
computed . 

2.  The  probe  separation,  Ay,  was  fixed  for  the  fold-over 
measurements  in  this  work.  It  would  be  interesting  to 
know  what  results  would  be  obtained  by  fixing  one  probe 
location  (preferably  at  the  mean  interface  location, 

Y)  and  varying  the  other  probe  location. 

3.  It  might  also  be  interesting  to  try  using  four  probes 
instead  of  two.  These  probes  can  be  located  in  pairs. 
One  pair  could  be  positioned  at  one  x/D  station  while 
the  other  pair  is  positioned  at  another  x/D  station. 

The  lateral  positions  should  be  such  that  two  lines 
joining  the  probes  are  parallel  and  form  an  angle  of 

45°  or  less  with  the  axis  of  the  jet.  This  angle  should 
be  measured  in  the  counterclockwise  direction.  However 

all  the  probes  should  be  perpendicular  to  the  flow. 

This  exercise  might  help  in  determining  the  shape  of 
the  interface,  and  possibly  the  width  of  the  fold-over 
region. 
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4.  The  effect  of  an  acoustic  disturbance  could  be  examined 
further  to  check  whether  the  percent  increase  in  fold- 
over  is  exactly  equal  to  the  percent  increase  in  the 
jet  widening  rate.  Since  this  work  considered  only 

one  Reynolds  number,  the  Reynolds  number  could  be  varied 
in  the  future. 

5.  It  was  pointed  out  earlier  that  the  jet  widening  rate 
decreased  with  increasing  Reynolds  number  at  frequencies 
of  measurable  sensitivity;  it  would  be  interesting  to 
know  whether  this  relationship  is  linear  or  not. 
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APPENDIX  A.  DIGITAL  PROCESSING  OF  DATA 

The  technique  used  to  digitize  the  tape  recorded 
analog  signal  will  be  presented  here.  A complete  block 
diagram  of  the  hybrid  computer  used  for  this  process  will 
also  be  given.  The  description  here  will  be  sufficiently 
detailed  so  as  to  initiate  future  experimenters  who  want 
to  use  the  system.  The  hybrid  computer  used  in  this  re- 
search is  located  at  the  Electrical  Engineering  building. 
This  computer  is  the  CDC  1700.  Before  the  actual  process 
is  described,  the  theory  underlining  discrete  representa- 
tion of  an  analog  signal  will  be  given. 

1 . Discrete  Representation  of  a Recorded  Analog  Signal 
Digitization  is  made  up  of  two  separate  and  distinct 
operations:  (a)  sampling  and  (b)  quantization.  Sampling 

is  the  process  of  defining  the  instantaneous  points  at 
which  the  data  are  to  be  observed.  Quantization  is  the 
conversion  of  data  values  at  the  sampling  points  into  num- 
erical form.  Sampling  for  digital  data  analysis  is  usually 
performed  at  equally  spaced  intervals  (see  Figure  A-1). 

The  sampling  interval  is  represented  by  h.  A problem 
arises  in  how  to  choose  h.  On  one  hand,  sampling  at  points 
which  are  too  close  together  will  yield  correlated  and 


r 
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highly  redundant  data,  and  thus,  unnecessarily  increase  the 
labor  and  cost  of  calculations.  On  the  other  hand, 
sampling  at  points  which  are  too  far  apart  will  lead  to 
confusion  between  the  low  and  high  frequency  components 
of  the  original  data.  This  latter  problem,  called  aliasing 
(see  Figure  A-2),  constitutes  a potential  source  of  error 
which  does  not  arise  in  direct  analog  data  p»'ocessing. 

This  will  be  discussed  further,  when  errors  associated 
with  digital  data  analysis  will  be  given. 

A continuous  random  record  x(t)  to  be  sampled  at 
points  At  = h apart,  can  be  seen  in  Figure  A-1.  The 
Nyquist  sampling  criterion  is  given  by  the  relation; 

h < 2 V (A-1  ) 

c 

where  F is  the  Nyquist  folding  frequency  (see  Bendat  and 
c 

Piersol  (1971)).  Two  practical  methods  exist  for  handling 
the  aliasing  problem.  The  first  is  to  choose  h sufficiently 
small  so  that  it  is  physically  unreasonable  for  data  to 
exist  above  the  associated  cutoff  frequency,  F^ . The  second 
method  is  to  filter  the  data  prior  to  sampling  so  that  in- 
formation above  a maximum  frequency  of  interest  is  no  longer 
contained  in  the  filtered  data.  Then  choosing  F^  equal  to 
the  highest  frequency  of  interest  will  give  accurate  results 
for  frequencies  below  F^ . This  second  method  is  recommended 
since  it  saves  on  computing  time  and  costs. 
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Let  represent  a set  of  data  values  found  at  the 

N points  where  n = 1,2,  ...,  N.  Then 

t^  = tQ  + nh  for  n = 1,  2,  N 

where  tQ  is  the  initial  point  chosen  and  N is  the  sample 
size.  This  sample  size  is  selected  on  the  basis  of  the 
desired  accuracy.  The  analog  signal  can  now  be  represented 
by  the  set: 

= X(tQ  + nh)  n = 1 , 2,  . . . , N (A-2) 

The  record  length,  T must  satisfy  the  equation 

T = Nh  (A-3) 

2 . Analog-Digital  Conversion  Technique 
The  description  of  the  technique  used  to  convert  an 
analog  signal  to  a digital  representation  is  presented  here. 
Care  is  taken  to  describe  the  steps  to  be  taken  starting 
from  the  laboratory,  where  the  data  is  tape  recorded,  to  the 
hybrid  computer  facilities,  and  finally,  to  the  computing 
center  at  the  Mathematical  Science  building.  Researchers 
at  Purdue  University  who  intend  to  digitize  their  data  and 
process  them  in  the  computer  will  find  this  section  extrem- 
ely valuable. 

In  the  last  section,  mention  was  made  about  the  method 
to  use  to  avoid  the  problem  of  aliasing  and  yet  to  be 
economical.  Using  that  method,  it  would  be  necessary  to 
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first  filter  the  signal  before  tape  recording.  In  this  j 

research  the  maximum  frequency  of  interest  was  found,  using  I 

I 

a spectral  analyzer,  to  be  below  2 KHz.  To  see  what  ef-  | 

feet  filtering  the  signal  had  on  the  digitization  process,  ' 

two  signals  were  recorded.  One  was  low  passed  at  2 KHz 
and  the  other  was  not.  No  difference  was  observed  after 
the  digitization  process.  This  was  because  the  maximum 
frequency  of  the  signal  as  obtained  by  the  spectral  analysis 
was  below  2 KHz.  It  was  not  necessary,  therefore,  to  fil-  l 

terthesignalbeforedigitizing.  i 

The  tape  recorder  used  was  an  FM  tape  recorder  model  | 

PI6104  (see  Figure  3-4  in  Chapter  III).  This  was  cali- 
brated before  any  recording  process.  It  should  be  pointed 
out  that  the  hybrid  computer  used  was  capable  of  quantizing 
voltages  in  the  range  -lOV  £ _<  lOV.  Any  signal  outside 

this  range  will  be  treated  as  -10  volts.  The  hybrid  com- 
puter assigns  a six  digit  number  to  the  full  scale  value. 

The  sign  (-)  is  included  in  the  number.  Specifically, 

-10  volts  = -32764.  However  if  the  signal  is  positive, 
the  sign  is  assumed.  If  the  nominal  voltages  of  the  signal 
were  significant,  they  can  be  obtained  by  graphing  the 
computer  numbers  against  voltages.  The  computer  number  for 
zero  voltage  is  zero.  It  is  necessary,  therefore,  to 
calibrate  the  tape  recorder  so  that  on  play  back  the  volt- 
ages fall  appreciably  within  -10  volts. 
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The  block  diagram  of  the  hybrid  computer  system  used 
in  the  digitization  process  is  shown  in  Figure  A-3.  A 
simplified  version  of  this  figure  was  shown  in  Chapter  III 
(Figure  3-5).  A description  of  the  patching  process  is  ; 

given  using  Figure  3-5. 

The  FM  tape  recorder  used  to  record  the  signal  should 
be  used  here.  The  two  channels  were  connected  to  the 
interface  rack  labeled  Trunk  070  and  071,  using  the  avail- 
able patch  cords.  The  Analog  trunk  lines  labeled  EAI  680 
P]  should  be  used.  This  is  internally  connected  to  the 
digital  computer.  On  the  Patch  Board,  ADC-0  was  con- 
nected to  Trunk  070  and  ADC-1  to  Trunk  071,  using  the 
appropriate  cords.  Also  on  the  Patch  Board,  a special 
logic  cord  was  used  to  connect  Sense-line-0  to  Pushbutton 

] 

^5.  The  patching  process  was  now  complete.  However,  the  i 

signal  could  be  viewed  on  the  oscilloscope  provided  for 
the  purpose,  since  one  of  the  oscilloscopes  is  internally 
connected  to  the  FAI  680  ^1  Analog  computer.  The  patching 
|j  could  be  checked  by  connecting  Trunk  070  or  071  to  Trunk 

l';  110.  The  signal  should  appear  automatically  on  the  scope 

1 if  the  scope  was  turned  on  and  the  tape  recorder  playing. 

The  hybrid  computer  consultants  should  be  consulted  in  case 

j of  problem. 

« 

j The  rest  of  the  process  involved  reading  a set  of 

I 


cards  and  manipulating  the  CDC  1700  console.  The  Teletype 
(1711)  was  used  to  give  command  to  the  rest  of  the 
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I 


equipment.  A step-by-step  procedure  on  manipulating  the 
CDC  1700  console  is  available  at  the  console  shelf. 
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The  program  used  to  digitize  the  analog  signal  was 
due  to  the  Hybrid  computer  consultants.  The  program  was 
capable  of  digitizing  up  to  4 channels  of  analog  data  using 
the  1538  in  conjunction  with  the  1571  (see  Figure  A-3).  The 
program  was  modified  to  accomplish  the  present  task.  Figure 
A-3  will  be  used  to  describe  the  logic  used  in  the  program. 

The  hybrid  system  was  capable  of  sampling  at  40  KHz. 

The  actual  sampling  rate  was  determined  by  the  number  of 
analog  channels,  NCtl,  and  the  number  of  5 microsecond 
intervals  per  sample,  NS.  For  example:  Assume  that  two 
analog  channels  of  data  will  be  digitized,  and  that  the 
maximum  frequency  of  interest  was  2 KHz.  The  sync  pulse 
interval  is  equal  to  the  h used  in  equation  (A-1),  So: 
h = sync  pulse  interval 

i .e. , h = NS  X 5 u s (A-4 ) 


■3 


Applying  equation  A-1: 


h 


2 X 10^ 


(A-5) 


Since  F^  was  already  determined  ' the  laboratory  to  be 
2 KHz,  NS  can  be  computed.  This  value,  50,  is  plugged  into 
the  program  and  the  sampling  rate  per  channel  will  be  1/h 
or  4000  samples  per  second.  Two  buffers  of  length  2500 
words  were  used  to  store  the  data  before  writing  on  the 


» 


I 

i 

I 

! 


i 

I 

i 


[i 

!■ 


magnetic  tape.  Since  the  two  analog  channels  were  digi- 
tized simultaneously,  a buffer  of  2500  words  will  contain 
1250  words  per  channel.  It  is  important  to  know  how  the 

data  is  arranged  in  the  user's  buffer  (the  magnetic  tape).  < 

Since  two  channels  were  digitized  using  ADC  #0  and  ADC  #1, 

the  output  will  read  ADC  n*!  before  ADC  #0.  However,  if  4 

channels  were  digitized  using  ADC  #0,  ADC  #1,  ADC  #2  and 

ADC  #3,  the  system  will  put  the  last  unit  digitized  first, 

i.e.,  ADC  #3  and  the  rest  in  series.  So  the  output  should 

be  interpreted  as:  ADC  n,  ACD  #0,  ADC  #1  and  ADC  #2.  Other 

subroutines  used  by  this  program  are  explained  in  the  hybrid 

program  documentation. 

A second  program  was  required  to  reformat  the  digi- 
tized data  for  CDC  6500  input.  This  program  can  also  be 
obtained  from  the  consultants. 

3 . Errors  Associated  with  Digital  Analysis 
In  section  one,  it  was  mentioned  that  one  of  the  prob- 
lems encountered  in  digital  data  analysis  was  that  due  to 
aliasing.  Figure  A-2  was  drawn  to  illustrate  this  problem. 

Assuming  the  sample  interval  was  chosen  to  be  h,  say,  the 
sampling  rate  will  be  1/h  samples  per  second.  However,  at 
least  two  samples  per  cycle  are  required  to  define  a fre- 
quency component  in  the  original  data.  Hence,  the  highest 
frequency  which  can  be  defined  by  sampling  at  a rate  of  1/h 
samples  per  second  is  l/2h  cps.  Frequencies  in  the  original 


data  above  l/2h  cps  will  be  folded  back  into  the  frequency 
range  from  0 to  l/2h  cps,  and  be  confused  with  data  in 
this  lower  range.  This  describes  how  aliasing  arises  (see 
Figure  A-2).  The  methods  of  handling  the  aliasing  problem 
were  described  in  section  one  of  this  Appendix.  Other 
errors  of  importance  are:  (a)  Aperture  error:  This  arises 
from  the  fact  that  the  data  sample  is  taken  over  a finite 
period  of  time  rather  than  instantaneously;  (b)  Jitter: 

This  arises  from  the  fact  that  the  time  interval  between 
samples  can  vary  slightly  in  some  random  manner;  and  (c) 
Nonlinearities:  This  is  due  to  many  sources  such  as  mis- 

alignment of  parts,  bit  dropouts,  quantization  spacing, 
zero  discontinuities,  etc. 

The  errors  a,  b,  c mentioned  above  cannot  be  corrected 
by  the  experimenter.  It  should  be  pointed  out  that  digi- 
tization process  merely  approximates  a continuous  analog 
signal  to  discrete  points.  Further  discussion  can  be  found 
in  Bendat  and  Piersol  (1971). 
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APPENDIX  B.  THE  FORTRAN  IV  PROGRAMS 

The  computer  programs  used  to  analyze  the  digitized 
analog  data  are  presented  and  described  here.  It  is 
obvious  that  after  a digitization  process,  a means  of 
evaluating  the  quality  of  the  data  is  essential.  Two  pro- 
grams were  used  for  this  purpose.  They  were  coded  SKIPS 
and  MSQ.  SKIPS  displays  any  portion  of  the  tape.  This  is 
accomplished  by  specifying  the  lines  at  which  to  start 
and  end.  Notice  that  in  the  control  card  that  requested 
the  tape,  the  parameter  C was  set  to  60.  This  was  because 
the  data  was  in  a binary  coded  format  which  used  6 bits  to 
represent  one  character.  There  were  ten  characters  per 
line  in  the  tape.  The  F = E is  used  to  describe  the 
type  of  format  to  be  used  (external).  A regular  size  mag- 
netic tape  can  store  up  to  90  buffers  of  words.  Each  buffer 
contains  exactly  2500  words.  Since  42960  points  were  an- 
alyzed for  each  location  (x/D,  y/b),  in  the  flow  field,  a 
little  over  17  buffers  will  be  required  to  hold  these  data 
points.  However,  since  two  channels  of  data  were  digitized 
simultaneously,  17  buffers  would  contain  S-1/2  buffers  per 
channel.  This  would  in  effect  mean  that  only  21480  data 
points  were  analyzed.  So  in  order  to  analyze  42960  words 


por  channel,  34.4  buffers  would  be  needed.  Since  a iiKuinetic 
t p e can  hold  up  to  90  t)  u f f c r s , it  can  be  seen  that  o ti  1 y 
four  locations  in  the  flow  field  can  be  analyzed  per  mag- 
netic tape.  A means  of  separating  the  data  foi-  each  pair 
of  location  atui  to  check  the  goodness  of  data  points  was 
developed.  Since  the  properties  of  a sine  or  triatigular 
wave  were  known,  either  of  these  signals  was  used.  Hence 
in  a magnetic  tape,  the  first  buffer  was  filled  with  either 
a s i tie  wave  or  a triangular  wave,  and  the  negative  full 
scale  values.  In  the  buffer  containing  TSOO  words,  12S0 
words  from  AOC  1 would  be  either  a sine  wave,  or  a tri- 
angular wave,  depending  on  which  signal  was  connected  to 
Trunk  line  071.  The  other  12S0  words  would  be  from  ADC  #0, 
depending  on  what  signal  was  connected  to  Trunk  line  070. 
Usually  nothing  was  connected  to  Trunk  line  070.  This 
would  have  the  effect  of  digitizing  the  negative  full  scale 
value  of  -32764.  The  next  44  buffers  contain  data  at  two 
locations  in  the  flow  field,  say  ( x / 0 = 35,  y ^ = 3.5  in 
(U.«9  cm),  y^  = y^  - Ay  - 3.5  - .125  = 3.375  in  (8.5725  cm)). 
Recall  that  the  probe  separation  was  Ay  - 0.125  in  (0.3175 
cm).  Probe  ^1  was  located  at  y^,  and  [irobe  f?  was  located 
at  y,,  , closer  to  the  Jet  center- line.  So  44  buffers  will 

L. 

contain  22  buffers  for  y^  and  22  buffers  for  However, 

0 n 1 y 4 2960  points  (about  17.2  b u f f e r s ) w t*  r e p r o c e s s e d per 
channel  as  was  already  noted.  The  other  test  signals  were 
c 0 ti  t a i n e d in  the  46 1.  h buffer,  and  finally  4 4 m o r e buffers 
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contained  the  other  pair  of  location  in  the  flow  field. 

These  second  test  signals  separate  the  first  set  of  data 
from  the  second. 

The  FORTRAN  IV  program,  SK1P3,  lists  any  portion  of 
the  magnetic  tape.  The  line  number  is  also  printed.  From 
the  printed  output,  ttie  data  were  examined  to  see  if  an 
overlap  of  data  existed.  This  is  accomplished  by  noting 
where  the  test  signals  started  and  stopped.  Figure  B-1 
is  a print  out  of  program  SKIP3. 

The  mean  square  of  the  voltage  signals  were  computed 
using  the  program  MSQ.  These  mean  square  values,  SI  and 
S2 , were  required  for  the  programs  that  constructed  the 
i n term i t tency  function,  I(t),  and  computed  the  fold-over 
events.  Due  to  the  limitation  in  the  core  memory  of  the 
CDC  6400-6500  dual  MACE  system,  (maximum  core  memory  150K), 
it  was  necessary  to  break  the  program  into  four  parts,  as 
given  by  NTIMES  = 4 in  the  program  MSQ  (see  Figure  B-2). 

This  program  also  located  the  maxima  for  each  set  of 
42960  points.  The  ratio  of  the  maxima  to  the  mean  square 
was  computed.  As  discussed  in  Chapter  III,  this  helped 
to  check  the  goodness  of  the  data  points.  If  the  maxima 
equaled  the  square  of  the  full  scale  values  (-32764)^,  it 
would  be  obvious  that  the  data  were  bad.  The  location  at 
which  the  maxima  occurred  was  also  printed.  With  this 
information,  data  points  surrounding  this  point  were  printed 
using  SKIP3.  Visually  checking  these  points  will  show  where 


the  data  became  bad. 
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After  checking  the  data  points  and  computing  the  mean 
square  values,  SI  and  S2,  the  next  step  involved  finding 
the  correct  threshold  level  to  use  in  order  to  construct 
the  "correct"  i ntermi ttency  function,  I(t).  This  exercise 
is  performed  using  the  program  THOLDIO.  The  ratio  of  the 
thresholo  level  to  the  mean  square  was  given  by  P = THl/Sl 
or  P = TH2/S2  (see  Figure  B-3).  This  program  varies  P at 
will.  For  each  P,  the  I(t)  function  is  constructed.  The 
hold  time  was  varied  by  the  equation: 

MM  = KK  - NT  (B-1 ) 

where  KK  denotes  the  numL^r  of  points  examined,  and  NT  de- 
notes the  number  of  points  that  govern  the  hold  time 
parameter.  If  NT  = 0,  the  hold  time  will  be  zero  and  if 
NT  = 40,  the  hold  time  will  be  computed  by: 

Th  = nrI^  ""  = 10  msec  (B-2) 

where  NRATE  is  the  sample  rate  which  was  determined  to  be 
4000  samples  /sec.  Notice  that  if  NT  = 20,  the  hold  time, 
will  be  equal  to  5 msec.  This  was  how  the  program 
varied  When  the  correct  was  obtained,  10  msec,  NT 

was  set  equal  to  40  for  the  rest  of  the  experiment. 

The  i ntermi ttency  function,  I(t),  was  formed  within 
the  "DO  LOOP": 

DO  20  N = 1,  MM  (B-3) 


(see  Figure  B-3).  Essentially,  the  program  compares  a given 


1 ‘'1  I I ■ I.  I'l' •llOl"'!' 1 , 'I'l'ii'-  I 1 L ‘ • t . I • '1. 

I .11. 

l-  t I > , •,  •ll  ll.l  > 

!|l"l|  I'  , liil.l  !*•.  U ll'.'- V'' ■ • I'-l  L ^ 

I'l  Ull..'-  I till  ill,-  ' 

I Hull  : i.;M-i.ii.i  1 11, • ' 

1 .1  I 1 1 1 1 . . M iiii  iH,-:  1 

•I  UK- 

n■■|ll  I nil  1lil.ll  fil  O'  1;IM|-Ur  IlIhU.''  h.fl  .-11  luGf, 

1 Inn  I.  nii|>'ii|  ' 

(HP  PF-iiiji-Yin  I IP;  n I lie:  l,k  n fii'.iiicnr  th' i « 

1H1H  PL  Hi',  null  rt-iii.L::t  . I'l.ipiu  i-i  'iiiirp:  ihpiihC' 
I L-ij.-i  tin  i itii:. 

SUi.i:i;i  I I'.i-fC  : :i  Ii-jti'  i'.  lb 

l I ll-ll  I.  iH  I ' 1 ij  ' > - 1 ■ 1 .1 1 1 1 ' • 1 U I'iJ  Cl  ■' 

F'^=-i:i.^ 

i.ij 

f-P+.5 

IP' HI  r.Ci?.  1 ''  GG  111  7'1 
PLUIHIi  3 

TM  HI  U(t;=300 

in  in  1 = 1 •Hi.  I HE 
LnLl-  FEnl'lKo' H Jl-lGH£>  l"k> 
iru  .LE.  O'  GU  TO  ipl 
1 0 C OHT 1 HUE 


CHtlii-E  THE  Pnil'-v  TO  111.  iPOCESSEl  UITH  [iE'nREr,i  VOLUE-?--  FV 
SPt  I.  lE'i'lHG  ISTHftT  mill  ILHP. 

HL  rH(  =E'l5E 
J1  = 0 
.ic  = 0 
t 1=0. 

KE'=0. 

HI  II1E  -S=4 

liij  , '■  fiElliiHi  -l"Hf  Hit  .1 
H I.  liPHI  .E'  . 1 ' CO  HI  40 
HI.  I;  . lEHli 
\ Ti.Pl--HI.  IHE*  l 
H H-.i  til.  IHL'HII  IHE 
til  - 1 . li-'i-r 
-il  Hfi 
1.1 1 11)  41 

4 0 1 . ri  .”  i '-HI  IHI"  1 1 

HI'  Hi  i HI.  H 11  1 1 lE 
Hi-i  |■|.;•■^ 
r.-  H'liii 

CdH  III  I Hi'  Vh'  ImH-  t 
U ; ' 

11' 

,''l  H I t1l'M.‘ 

Pi  i.:  ' • ■ i ■ H ' • H I • I 0 ' 

'■11  I ■ H I • iri.p 
I I 1 I I 1 

I ■:  i I I ■ I'l  p,  I I ■ It  ' • 

,'  I 1 i H , ■■  n -p 
I I -1  I ' 1 

14  ' . I . ■ ) I H.H  ' 1 ' H ' • ■ 

Figure  B-3.  Program:  THOLDIO. 


iiMiiriiif 

3 ■:  I 1 .1  . . 1 1 :i  1 1 ■ 'I 

■ ■1  I.'  • 

; I-  , f , ■ ■ 

'S  . ■ ’ 1 I ■ ■ ■ • •;  li!  t-'i.l.’i  ' ■> 

•r.'.  I nril  ill'.if 

! ■ '.I  t,  ■!  M.’ 

■ ■ ■ 

1 H I t ' ■ 1 

nr=4i:i 

i. -MT 

L=tKint:vrtM 

I ii*4i  i i.-4  i.tOO 

T “Fl  [jmT  ‘v  1 1 ' f 1 PAT Pf-'H T ' 

l^  Uli,  ui.rtT.til  . 1'  iVI  Uj  4 3 
mp;  I c • 1 " 

s tuKni-)r<-  • ‘-V  . I [ii  I .)L  s.'.urLr.,  i-i;''  ; i ■.  'ifLU-iG-f 

1 iF'i'.:' ^Tiir.iL  I ctf.  j’-rc--..;  .••■..-ff.  < ■'.■..Ps:  :i  -• 

■iFlS.  .("JuKLIc  i-.. . lO.- • ■'THPE:  tiPi.li  !:■  hK 

43  ri;l  30  H=i.rin 

:f*ri.riF.  1 ' c.n  to  pi 

UWiaiJ.uK.THl''  niJ  TD  33 
DQ  on  .1.1=1,. to 

30  IF' VI ■ lifjj'.i’it:. Trii Gi]  rn  33 
Oil  I U 33 

33  = 

.U  = .ll4l 

iiil  iO  34 

33  ‘v'i'rF'  = 0. 

34  1F<lV3>:M.>.GF..TH3>  Gil  TO  ST. 
lill  01  Jj=1.40 

1 4 \ V.- ' 11+- JJ ' , OF , 1 113 ) oil  in  cT 

G;  i i li  30 

3T  v3<.‘F--l. 

JS-J.'.'-t  1 
Gil  Til  30 
,3>-  V.3  <0  . = 0 . 

GJ  ill  30 

31  ..  r • 1,1  .■■7 

IF  • (M-l  ' . OF.  0.  > : if  33 

'.1  ■ it  ■ - 0. 

on  111 

■>:  nil  I .1.1-1 ,41: 

I'c  14  ' '1  i’4.i.M,!.;f  . TH I ' Oil  ID  .JF 
GU  111  ;0 


.1  ".lit  1 
I'.J  III 

50  Vl'ti'-O. 

4 r 4 : : • . 

I-'  ti  ' r.  .GL.  • !.'■  lU  ■ ; 
ll  ' . . 0.  ' rij  , > 

V. ' O ■ • I.' . 
l.-l  1.  . .1 

53  liJ  4-  ■ . .'•t<  111 
• 5 ll'  >.  I .1.'  .'F.Tn.'i  I'll  111  ' I 
I p 1(J 

51  V 1 . 

1,  I.  ' I 
I'll  i.i  0 
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Figure  B-3.  Continued. 
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point,  V1(N)  to  the  threshold  level,  Ttll.  If 

VI (N)  > THl , VI  = 1 (B-4) 

Otherwise  V1(N)  will  not  be  allowed  to  change  level  until 
it  passes  the  hold  time  criterion.  This  criterion  checks 
40  future  points  to  see  if  any  of  them  is  greater  than  THl. 
If  all  of  the  40  points  are  each  less  than  THl,  V1(N)  will 
be  set  equal  to  zero.  Recall  that  the  40  points  were  tan- 
tamount to  a hold  time  of  10  msec.  Notice  from  the  program 
that  before  the  hold  time  criterion  is  applied  to  the 
point,  V1(N),  it  first  checks  the  value  assigned  to  Vl(N-l). 
Obviously,  if  Vl(N-l)  = 0,  V1(N)  would  automatically  be  set 
equal  to  zero.  This  is  because  in  order  to  set  Vl(N-l) 

= 0,  V1(N)  was  found  to  be  less  than  THl.  The  hold  time 
criterion  is  applied  to  Vl(N)  only  when  Vl(N-l)  = 1.  The 
obvious  question  would  be:  What  happens  when  N = 1?;  that 
is,  when  the  first  point  is  examined  since  the  computer 
"blows  up"  at  V1(0).  The  program  accounts  for  this  by 
checking  if  N = 1.  This  time  it  compares  VI (N),  which  is 
Vl(l),  to  THl.  If  Vl(l)  > THl,  Vl(l)  = 1.  Otherwise,  the 
hold  time  criterion  will  be  applied. 

The  program,  THOLDIO,  computes  > and  for  each  value 
of  P.  As  was  discussed  in  Chapter  III,  curves  of  'i  and 
vs.  P are  drawn  to  determine  the  range  of  P where  F^  re- 
mains relatively  constant  and  begitis  to  approach  an 
asymptotic  value.  From  these  curve,  the  "correct"  P is 


chosen . 
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As  soon  as  the  "correct"  P is  selected  (actually  PI 
and  P2  for  channels  1 and  2,  respectively)  the  values  are 
plugged  into  the  fold-over  program,  FOLD.  The  logic  in 
choosing  the  type  of  fold-over  event  was  discussed  in 
Chapter  III  and  will  not  be  repeated  here.  A listing  of 
program,  FOLD  is  shown  in  Figure  B-4.  Tables  B-1  and  B-2 
show  the  values  used  to  plot  y and  F^  vs.  y/b  for  the 
undisturbed  and  acoustically  excited  cases,  respectively. 
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Figure  B-4.  Program;  FOLD. 
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G£  II ' ,1.1  .or,.  THl  ' GO  TO  £■" 

Go  TO  50 
£7 
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Figure  B-4.  Continued. 
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i -■  Iji  Jl-*.- ; 'f  I ILD4  + 1 

tn  ,ri  1=0 
nr, I t £-0 
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Figure  B-4.  Continued. 
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Figure  B-4.  Continued 


Table  3-1.  Values  for  7 and  F vs.  y/b:  Undisturbed  Case. 
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Table  B-2,  Values  for  y and  F vs.  y/b:  Acoustically  Excited  Jet. 
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Table  B-3.  Percent  Fold-Over:  Undisturbed  Case. 


x/D 

y/b 

Y 

Fro?it 

Bak 

Averlge 

0.913 

0.978 

16.67 

27.78 

22.22 

1.52 

0.617 

40.35 

33.14 

36.73 

35 

1.83 

0.294 

29.56 

39.38 

34.48 

2.13 

0.136 

43.30 

39.80 

41.54 

2.74 

0.015 

63.16 

78.95 

71.05 

0.714 

0.973 

22.73 

36.36 

29.55  r 

1.19 

0.851 

30.56 

16.67 

23.61  S 

45 

1.66 

0.446 

28.78 
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34.77  U 

2.14 

0.127 

64.18 

62.12 

63.16  i 

2.62 

0.015 

80.85 

81.25 

81.05  f 

0.98 

0.933 

30.61 

40.82 

35.71 

1.37 

0.763 

45.51 

46.84 

46.18  1 

55 

1.17 

0.850 

20.83 

30.14 

25.52  i 

1.76 

0.306 

55.83 

54.60 

55.21  ' 

2.15 

0.098 

67.62 

71.43 

69.52 
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Table  B-4. 


Percent  Fold-Over: 
spl  = 105  dB, 


Acoustically  Excited  Case:  F = 1149  Hz, 


x/D 

y/b 

Y 

FroAt 

Ba^k 

Average 

1.07 

0.867 

29.33 

30.26 

29.80 

35 

1.61 

0.465 

38.27 

40.99 

39.63 

1.88 

0.235 

43.75 

49.56 

46.67 

1.06 

0.869 

39.44 

31.94 

35.66 

45 

1.48 

0.522 

54.93 

47.55 

51.22 

1.90 

0.155 

48.15 

66.67 

57.41 

1.05 

0.798 

31.58 

33.33 

32.46 

1.22 

0.733 

34.65 

36.89 

35.78 

1.57 

0.329 

57.14 

55.63 

56.38 

55 
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APPENDIX  C.  THE  DIFFERENTIATOR  CIRCUIT 


The  differentiator  circuit  used  to  differentiate  the 
hot-wire  voltages  is  described  here.  When  the  circuit  was 
constructed,  a known  signal,  sine  wave,  was  passed  through 
it  to  see  if  a cosine  wave  was  constructed.  When  a satis- 
factory combination  of  resistors  and  capacitors  was  found, 
the  values  were  held  constant  for  the  rest  of  the  experi- 
ment. 

The  differentiator  circuit  used  is  shown  in  Figure 

C-1 . 

The  input  impedance,  Z^,  is  given  by  the  equation: 


Z,(u)) 


(C-1  ) 


Where  Rj  and  c are  the  resistor  and  capacitor,  to  is  fre- 
quency, and  j is  a complex  number  defined  by  j = /H".  The 
feedback  impedance  is 


Zf(u)  = R2  (C-2) 

Using  the  Laplace  variable,  s,  equations  (C-1)  and  (C-2) 
can  be  written  as 


(C-3) 
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and 


Z^(s)  = R2 


(C-4) 


Using  this  notation,  the  Transfer  Function  for  the  circuit 
is  given  by 

e (5)  Zf(s) 

* irnr  ■ ■ zrrrr 


* cs 


(82^)5 
{R^C%  + 1 


(C-5) 


Letting  jw  = s , equation  C-5  becomes 

(R2C)  (jw) 

TF(u.)  = - (ju.)  “a 


(C-6) 


The  combination  R2C  is  known  as  the  constant  element,  j 
is  known  as  the  differentiator  element.  The  break  frequency, 
or  cutoff  frequency,  is  defined  by: 


1 

‘*>b  " rTC 


(C-7) 


Taking  the  magnitude  of  TF , we  get 
|TF(w)|  = 


R2C  w 


+ 1 

‘^b 


^/2 


(C-8) 


for  (i)  = lou.  we  get 


^^b 

/r 


1TF(U))| 
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For  u)  <<  we  get 

|TF(u))|  = R^c  u, 

For  u)  >>  u).  : 

0 

I TF  (u)  I = R^c 

A sketch  of  lTF(u))|,  showing  the  approximate  differentiator 

range  of  application,  can  be  seen  in  Figure  C-2.  It  can  be 

seen  that  for  frequencies  less  than  w,  the  circuit  will 

b 

act  as  a good  differentiator  circuit.  The  value  of  to,. 

b 

recommended  is: 

10^  ~ 10  X highest  frequency  of  interest. 

As  was  previously  mentioned,  this  highest  frequency  was 
2 KHz.  Hence,  ^ = 20  KHz.  The  values  of  R^,  R2  and 

c used  were:  R^  = 7.5  k,  R^  = 750  k,  and  c = 0.001  pf. 


Figure  C-2.  Sketch  of  |TF|  vs.  y/b 


